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Context. Measuring the mass of galaxy clusters is a key issue in cosmology. Among the methods employed to achieve this goal, the techniques 
based on lensing and X-ray analyses are perhaps the most widely used. However, the comparison between these mass estimates is often difficult 



Q and, in several clusters, the results apparently are inconsistent. 



Aims. We aim at investigating potential biases in lensing and X-ray methods to measure the cluster mass profiles. 



CZ) Methods. We do so by performing realistic simulations of lensing and X-ray observations that are subsequently analyzed using observational 
, ^ , techniques. The resulting mass estimates are compared among them and with the input models. Three clusters obtained from state-of-the-art 

hydrodynamical simulations, each of which has been projected along three independent lines-of-sight, are used for this analysis, 
r- H Results. We find that strong lensing models can be trusted over a limited region around the cluster core. Extrapolating the strong lensing mass 

models to outside the Einstein ring can lead to significant biases in the mass estimates , if the BCG is not modeled properly for example. Weak 
C^) lensing mass measurements can be largely affected by substructures, depending on the method implemented to convert the shear into a mass 
\l estimate. Using non-parametric methods which combine weak and strong lensing data, the projected masses within R 2 m can be constrained 

with a precision of ~ 10%. De-projection of lensing masses increases the scatter around the true masses by more than a factor of two due to 
' I cluster triaxiality. X-ray mass measurements have much smaller scatter (about a factor of two smaller than the lensing masses) but they are 
^SJ generally biased low by 5 - 20%. This bias is entirely ascribable to bulk motions in the gas of our simulated clusters. Using the lensing and the 
t-H X-ray masses as proxies for the true and the hydrostatic equilibrium masses of the simulated clusters and by averaging over the cluster sample 

we are able to measure the lack of hydrostatic equilibrium in the systems we have investigated. 

Conclusions. Although the comparison between lensing and X-ray masses may be difficult in individual systems due to triaxiality and sub- 
structures, using a large number of clusters with both lensing and X-ray observations may lead to important information about their gas physics 
• and allow to use lensing masses to calibrate the X-ray scaling relations. 

X 

^ Key words. Galaxies: clusters: general; X-ray: galaxies: clusters; Gravitational lensing 



1. Introduction 

Galaxy clusters are highly important test sites for cosmology. 
They are the most massive gravitationally bound structures in 
the universe and, in the framework of the hierarchical structure 
formation scenario, they are also the youngest systems formed 
to date. For this reason, the interplay between baryons and 
dark-matter, and its effects on the cluster internal structure, is 
less important in these than in smaller and older objects. Thus, 
they are ideal systems for testing the predictions of the cold- 
dark-matter paradigm on the internal structure of dark-matter 
halos (Yoshida et al. 2000; Clowe et al. 2004; Markevitch et al. 



2004; Sand et al. 2008). Moreover, their mass function is highly 
sensitive to cosmology, since its evolution traces with exponen- 
tial magnification the growth of the linear density perturbations 
(Press & Schechter 1974; Sheth & Tormen 2002; Jenkins et al. 
2001; Warren et al. 2006). 

For these reasons galaxy clusters are being used to con- 
strain the cosmological parameters. Cluster masses are used 
to measure the time evolution of the mass function, which is 
compared to the theoretical predictions to constrain the contri- 
bution of the components of the universe to the overall density 
parameter, the equation of state of dark energy and the normal- 
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ization of the power spectrum of the initial density fluctuations 
(e.g Vikhlinin et al. 2009; Mantz et al. 2009a,b). The measure- 
ments can be compared with those obtained from the observa- 
tion of the universe on larger scales, and combined with CMB 
(see e.g. Komatsu et al. 2009), Supernovae la (e.g. Riess et al. 
1998, 2004; Perlmutter et al. 1999), and Baryonic- Acoustic- 
Oscillations observations (e.g. Eisenstein et al. 2005; Percival 
et al. 2007), in order to put tighter limits to the values of the cos- 
mological parameters. A different approach consists of measur- 
ing the concentration-mass relation of clusters and of compar- 
ing it with the theoretical predictions (e.g. Pointecouteau et al. 
2005; Vikhlinin et al. 2005; Buote et al. 2007; Schmidt & Allen 
2007; Comerford & Natarajan 2007). For example, Buote et al. 
(2007) find that the concentration-mass relation measured with 
a sample of 39 clusters is good agreement with the expectations 
in the framework of the concordance model. 

Allen et al. (2008) use the cluster gas fraction to constrain 
the time evolution of the dark energy component of the uni- 
verse. The gas fraction measured within a given angular radius 
is proportional to the distance of the cluster to the power 1.5. 
Thus, the evolution of the gas fraction can be used to measure 
the cosmic acceleration (Allen et al. 2004). Combining the re- 
sults from the f gas technique with CMB and SNIa data sets, 
they find that the time evolution of the dark energy equation 
of state is compatible with the cosmological constant paradigm 
(see also Ettori et al. 2003). 

These techniques rely on scaling relations which link the 
mass to X-ray observables, such as temperature, pressure and 
luminosity of the X-ray emitting intra-cluster gas. The scaling 
relations are predicted to have limited scatter in mass. For ex- 
ample, simulations suggest that the mass-7x and the mass-Tx 
relations have just 15% and 8% scatter in mass at fixed Tx and 
Y x (e.g. Evrard et al. 1996; Kravtsov et al. 2006). Despite these 
encouraging predictions, it is clearly essential to measure and 
calibrate scaling relations empirically. To this purpose, we need 
to accurately measure the mass profiles of an as large as possi- 
ble sample of galaxy clusters. 

There are several methods to derive the masses of clusters. 
Two widely used approaches are based on X-ray and lensing 
observations. X-ray observations allow to derive the cluster 
mass profiles assuming that these systems are spherically sym- 
metric and that the emitting gas is in hydrostatic equilibrium 
(e.g. Henriksen & Mushotzky 1986; Sarazin 1988; Ettori et al. 
2002). This method has the advantage that, since the X-ray 
emissivity is proportional to the square of the electron density, 
it is weakly sensitive to projection effects due to masses along 
the line of sight to the clusters. However, it is still not well es- 
tablished how safely the hydrostatic equilibrium approximation 
can be made. 

Being the largest mass concentrations in the universe, 
galaxy clusters are the most efficient gravitational lenses on 
the sky. Their matter distorts background-galaxy images with 
an intensity which increases from the outskirts to the inner re- 
gions. Strong distortions, leading to the formation of "gravi- 
tational arcs" and/or to the formation of systems of multiple 
images of the same source, occur in the cores of some mas- 
sive galaxy clusters. Weak distortions, which can be measured 
only statistically, are impressed on the shape of distant galax- 



ies which lay on the sky at large angular distances from the 
cluster centers (e.g. Bartelmann & Schneider 2001). Both these 
lensing regimes can be used to map the mass distribution in 
galaxy clusters. Determining the masses and the density pro- 
files using lensing offers several advantages compared to X- 
ray observations. First, lensing directly probes the cluster total 
mass, including the dark matter component, without need to 
make strong assumptions on the equilibrium state of the lens. 
Second, mass profiles can be measured over a wide range of 
scales, from < 100 kpc out to the virial radius. The biggest 
disadvantage is that lensing measures the projected mass in- 
stead of the three-dimensional mass. It is much more sensitive 
than X-ray methods to projection effects, such triaxiality and 
additional concentrations of mass along the line of sight. Given 
the pros and cons of each method, we can conclude that lens- 
ing and X-ray are complementary to each other in many ways. 
In particular, the comparison of these two mass estimates can 
greatly help to improve the accuracy of the measurements and 
to understand the systematic errors. 

The picture arising from the comparison of lensing and 
X-ray mass estimates is puzzling. The two estimates are of- 
ten discrepant, which implies that we may be missing impor- 
tant ingredients for fully understanding the properties of galaxy 
clusters. A systematic discrepancy has been revealed in the 
sense that masses derived from strong lensing are typically 
larger by a factor 2-3 than masses derived from weak lensing 
and from the X-ray emitting intra-cluster-medium (ICM) (e.g. 
Wu & Mao 1996; Ota et al. 2004; Smith et al. 2005; Riemer- 
S0rensen et al. 2009). The comparison between weak lensing 
and X-ray mass estimates is less problematic (see e.g. Ettori & 
Lombardi 2003; Allen et al. 2001) but still the lensing masses 
are on average larger than the X-ray masses by ~ 15 - 20% 
(e.g. Hoekstra 2007; Zhang et al. 2008). Explaining this dis- 
crepancy was found to be difficult. It was proposed that X-ray 
masses could be systematically biased because of mergers or 
bulk motions of the gas that alter the state of hydrostatic equi- 
librium (Bartelmann & Steinmetz 1996; Allen 1998; Mahdavi 
et al. 2008; Piffaretti & Valdarnini 2008; Ameglio et al. 2009; 
Lau et al. 2009). As stated above, lensing masses could also 
be affected by significant uncertainties (e.g. Bartelmann 1995; 
Meneghetti et al. 2007). 

In this paper, we study the systematic effects in mass mea- 
surements done with standard lensing and X-ray techniques. 
Our work extends some previous works on the systematics on 
X-ray based mass measurements (see e.g. Rasia et al. 2006; 
Nagai et al. 2007), which used synthetic X-ray observations 
of clusters obtained from hydrodynamical simulations to re- 
cover their mass distribution under the assumption of hydro- 
static equilibrium. Here, we simulate both optical and X-ray 
observations of the same simulated clusters, in order to be able 
to compare the mass estimates obtained from both kind of ob- 
servations. We use three clusters obtained from hydrodynam- 
ical simulations and study them along three independent lines- 
of-sight. 

The paper is organized as follows. In Sect. 2 we describe 
the sample of numerically simulated galaxy clusters, as well as 
the techniques used to simulate lensing and X-ray observations. 
In Sect. 3, we introduce the methods through which the mass 
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estimates are derived from the mock data. In Sect. 4 we show 
the results of the analyses. Finally, in Sect. 5 we summarize and 
discuss the major findings of our study. 



2. Simulations 

In this Sect, we describe the numerical methods used 
in the simulations. We start by considering the N- 
body/hydrodynamical simulations from which the cluster mod- 
els are obtained. Then, we explain the lensing and the X-ray 
simulation pipelines. 



2.1. N-body/SPH simulations 

The clusters used in this work are gl, g5l, and g72 from the 
sample of numerical hydrodynamical simulations presented by 
Saro et al. (2006). These objects have been already used in sev- 
eral other studies (Dolag et al. 2005; Puchwein et al. 2005; 
Meneghetti et al. 2007, 2008; Rasia et al. 2006, 2008). Thus, 
we refer the reader to these papers for more details. They are 
extracted from a parent simulation of only dark matter (Yoshida 
et al. 2001) with a box size of 479/T 1 Mpc of a flat ACDM 
model with Q. m = 0.3 for the present matter density parame- 
ter, h = 0.7 for the Hubble constant in units of 100 km 
Mpc -1 , cr 8 = 0.9 for the rms fluctuation within a top-hat sphere 
of 8/z _1 Mpc radius, and Q;, = 0.04 for the baryon density pa- 
rameter. Each of them is subsequently re-simulated at higher 
mass and spatial resolution using the Zoomed initial condition 
method (Tormen et al. 1997). The re-simulations are carried 
out with GADGET-2 1 (Springel 2005). In the code we include 
(i) a description of several physical processes in the ICM (see 
Saro et al. 2006), (ii) a numerical scheme to suppress artificial 
viscosity far from the shock regions (see Dolag et al. 2005), and 
(iii) a treatment of chemical enrichment from both SNIa , SNII, 
as well as from low and intermediate mass stars (Tornatore 
et al. 2004, 2007). Our simulations assume the power-law 
shape for initial stellar mass function, as proposed by Salpeter 
(1955), and galactic ejects with a speed of 500 km s _1 . They 
start with a gravitational softening length fixed at e = 30/r'kpc 
comoving (Plummer-equivalent) and switch to a physical soft- 
ening length of e = 5/z~'kpc at z = 5. The final masses of the 
DM and gas particles are set to m DM = 1.13 x 10 9 /z _1 M o and 
wigas = 1.7 X 10 8 /i _1 M Q , respectively. 

During the simulation 92 time slices are saved from redshift 
60 to 0. These are equidistant in time. For the current work, 
we use the snapshots corresponding to redshift z\ = 0.297 for 
gl and g72 and to redshift zi = 0.2335 for g5l. These red- 
shifts are optimal for both the X-ray and the lensing analyses, 
since 1) the X-ray surface brightness of these objects is high, 
and 2) the strength of the lensing signal for sources at z s ~ 1 , 
which is approximately the median redshift of the sources in 
our simulations, is maximal at z\ ~ 0.3. All clusters are mas- 
sive, with masses M200 ranging between ~ 7 x 10 14 hr l M Q and 
~ 1.2x 10 15 hr l M Q (see Table 1). M200 is the mass enclosed by 

1 http://www.MPA-Garching.MPG.DE/gadget/ 



the radius r2oo, i.e. the radius within which the average density 
is 200 times the critical density, 

3// 2 (z) 



Per = 



8ttG 



(1) 



where H(z) is the Hubble parameter. 

For each cluster, we select a cube of 20/z~ 1 Mpc side length 
(comoving), centered on the most bound dark matter particle. 
This is sufficiently large to cover a wide field of view, needed 
for weak lensing simulations. All the matter contained into this 
box is projected along three orthogonal lines of sight, in order 
to produce three lens planes for each cluster. The surface den- 
sity maps corresponding to the three projections of the clusters 
gl, g5l, and g72 are shown in Fig. 1. The corresponding mass 
profiles are displayed in Fig. 2. Note that the sample of clusters 
considered here comprises object with different morphologies, 
triaxialities, and levels of substructures, although in the X-ray 
they all appear quite relaxed (see Fig. 4). While the cluster gl 
has the most regular morphology, the cluster g72 has a mas- 
sive companion (~ 10 14 h~ l M Q ) located at ~ 2.5 /r'Mpc from 
the cluster center. In the projection along the z-axis, this sec- 
ondary clump is much closer to the main halo (~ 300 /z _1 kpc). 
Also, the cluster g5 1 has some substructures within the inner 
500 /r'kpc but their masses are smaller (<5x 10 12 h~ l M Q ). 
A massive clump of mass 5 x 10 13 lr l M Q orbits at a distance 
of 3 /z _1 Mpc from the center. 

As demonstrated by the differences between the 2D-mass 
profiles, all three clusters are triaxial. Their shape is prolate, 
with the major axis oriented nearly along the z-axis of the sim- 
ulation box in the cases of g5l and g!2, and nearly along the 
x-axis in the case of gl. The axis ratios, measured at rxoo by 
calculating the cluster inertia ellipsoids, are listed in Table 1, 
together with some other relevant properties. We also report 
there the angles between the main axes of the clusters and the 
x-, y—, and z-axes of the simulation boxes. 

The clusters are well described by Navarro-Frenk-White 
(NFW) density profiles (Navarro et al. 1997), whose functional 
form is given by 

Ps 



pNFw(r) 



(2) 



r/r s (l +r/r s ) 2 

where p s and r s are the characteristic density and the scale ra- 
dius, respectively. The characteristic density is often written in 
terms of the concentration parameter, C200 = '"20o/'"s, as 



200 



"200 



Ps 3 Pci [ln(l + c 200 ) - caoo/(l + c 200 )] 



(3) 



The last two columns of Table 1 report the best-fit concentra- 
tions and scale-radii obtained by fitting the 3D dark-matter den- 
sity profiles of the clusters with the formula in Eq. 2. The radial 
fits are done between 10/z~'kpc and r2oo- 

In the remainder of the paper, we will refer to the projec- 
tions of cluster gN along the simulation axis X using the abbre- 
viation gN-X . 

2.2. Lensing simulations 

In this section we describe the ray-tracing algorithms used to 
derive the lensing distortion fields of our simulated clusters. 
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Fig. 1. Surface density maps of the three projections of the clusters gl (top panels), g5l (middle panels) and g!2 (bottom panels) used in this 
work. The size of each map is 10 /r'Mpc comoving. Such a scale corresponds to ~ 51.85 arcmin for g51 and to ~ 41.47 arcmin for gl and gll. 
In each panel, the inner, middle, and outer boxes indicate the fields-of-view used for the HST, Chandra, and SUBARU simulated observations, 
respectively (see text for more details). 



For a summary of the lensing definitions used through the rest 
of the paper, we refer the reader to Appendix A. 

In the following, we make use of the standard thin-screen 
approximation, i.e. we assume that the deflections occur on a 
plane perpendicular to the line-of-sight, passing through the 
cluster center. This is justified since the distances between the 
observer and the lens and between the lens and the background 
sources are much larger than the sizes of the clusters. 

The particles projected on each lens plane are used to cal- 
culate the deflection angles of bundles of light rays. The light 



rays are traced from the observer position towards the back- 
ground sources through two regular grids with different spa- 
tial resolutions. The inner 1.5 x 1.5 /T 2 Mpc 2 region around the 
cluster center is sampled with 2048 x 2048 light rays. This 
guaranties sufficient spatial resolution for reproducing accu- 
rately the positions of multiple images in the strong lensing 
(SL) regime (Meneghetti et al. 2007). For the weak lensing 
(WL) regime, we need to sample a much wider area, while the 
spatial resolution is less important. Thus, we cover the whole 
lens plane with a grid of 4096 x 4096 light rays. The deflec- 
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Table 1. Main properties of the simulated clusters used in this work. Column 1: cluster name; Column 2: redshift; Column 3: rioo; Column 4: 
M 2 (io\ Columns 5-6: principal axes ratios: b/a, c/a, where a> b > c; Columns 7-9: angles between the main principal axis and the x-, y—, and 
z-axes of the simulation box; Column 10: best-fit 3D-concentration; Column 11: best-fit 3D-scale radius 
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Fig. 2. Mass profiles of the clusters gl, g51, andg72. The solid black and red lines indicate the total and DM only 3D-mass profiles, respectively. 
The total 2D-mass profiles corresponding to the x, y, and z projections of each cluster are given by the dotted, dashed, and dash-dotted lines. 

tion angles are computed using a tree-based code, which works 1 " " A 



P V = uZ m « r|2 ^' (6) 



u=\ 



as follows. First, it ranks the particles based on their distances 'i ~ ~m 
from the light ray positions, building a Barnes & Hut oct-tree 

in two dimensions. The contributions to the deflection angles Assuming a Plummer softening to avoid that the deflection an- 

from nearby and distant particles are calculated separately us- gles diverge, the Fu{R') functions are defined as: 
ing direct summation or higher-order Taylor expansions of the 

deflection potential around the light-ray positions. Precisely, F\(R') = (7) 

given a light ray at position R in physical units, which corre- + s ) 

sponds to an angular position = R/D;, the contribution to its f 2 (^') = ~^ (g) 

deflection angle by a system of mass elements, m a at positions (R' 2 + s 2 ) 2 



R a (a = 1,2, ...,N - l,N), with center of mass Rcm, and with 
t\ « \h 

4GM 



R - R C m\ •« l*a - Rcm\ for all the mass elements a, is ^ (R'2 + s 2)3 



a t (R) 



1 ' - ~ 2 

F^R'W + F 2 (R')Q^ + -F 3 (R')(R' k Q k „R'"W J Fa{R } " {R 2 + s 2 ) 2 ' (10) 



\f 4 (R')p u 



2 

Nearby particles are treated as point lenses and Eq. 4 reduces 
R) (4) to 



4GM 

where M is the total mass of the system, /?'=/?- R CM , 6 lJ is »,(/?) = — —R^F^R ) . (11) 
the Rroneker function, and the tensors P and Q are defined as: c 



£1 = ^ 



The fraction of of particles which are evaluated with Eqs. 4 
q _ _ V 1 m R' a R' a (5) or 1 1 is set by the Barnes-Hut opening criterion, Obh (see e.g 

1^1 " ' ' ' Springel 2005), which we fix at 8 BH = 0.4. As shown by Aubert 
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et al. (2007), the optimal softening length s depends on the 
resolution of the simulation. We performed several tests to de- 
termine which values to use. Doing ray-tracing through NFW 
halos sampled with a similar number of particles as our simu- 
lated clusters, we verified that a softening scale of 5 /r'kpc is 
appropriate for reliably reproducing the deflection angle field 
of the input models over the range of scales relevant for both 
strong and weak lensing. 

Having obtained the deflection angle maps, we apply the 
cluster distortion fields to the images of a large number of 
background galaxies. While doing so, we simulate optical ob- 
servations of each cluster under the different projections. For 
this purpose, we use the code described in Meneghetti et al. 
(2008) (quoted as SkyLens hereafter), which has been recently 
further developed. In short, the code uses a set of real galax- 
ies decomposed into shapelets (Refregier 2003) to model the 
source morphologies on a synthetic sky. In the current version 
of the simulator, the shapelet database contains ~ 3000 galax- 
ies in the z-band from the GOODS/ACS archive (Giavalisco 
et al. 2004) and ~ 10000 galaxies in the B, V, i,z bands from 
the Hubble-Ultra-Deep-Field (HUDF) archive (Beckwith et al. 
2006) . Most galaxies have spectral classifications and pho- 
tometric redshifts available (Bemtez 2000; Coe et al. 2006), 
which are used to generate a population of sources whose lumi- 
nosity and redshift distributions resemble those of the HUDF. 

SkyLens allows to mimic observations with a variety of 
telescopes, both from space and from the ground. In this work 
we simulate wide field observations, on which we carry out a 
weak lensing analysis, using the SUBARU Suprime-Cam. We 
simulate Hubble-Space-Telescope observations of the cluster 
central regions using the Advanced Camera for Surveys (ACS). 
All simulations include realistic background and instrumental 
noise. As an example, in Fig. 3 we show two color-composite 
images of one simulated galaxy cluster (although in the rest of 
the paper we will use simulated observation in a single band). 
The left and the right panels show the results of simulated ob- 
servations with SUBARU and with HST, respectively. Note that 
the galaxy colors are realistically reproduced by adopting 22 
SEDs to model the background galaxies, following the spectral 
classifications published by Coe et al. (2006). 

In the HST image, several blueish arcs and arclets are vis- 
ible behind the cluster galaxies. These are originated by back- 
ground spiral and irregular galaxies strongly lensed by the fore- 
ground cluster. 

In our analysis, we use simulated observations with the fol- 
lowing characteristics. For the SUBARU simulations, we as- 
sume an exposure time of 6000s in the / band, with a seeing 
of 0.6". The PSF is assumed to be isotropic and modeled us- 
ing a two-dimensional Gaussian. For the HST simulations, we 
assume an exposure time of 7500s with the F775W filter. The 
fields-of-view adopted for the HST and SUBARU simulations 
are overlaid to the surface density maps in Fig. 1 (blue inner- 
and white outer-boxes, respectively). Note that these fields of 
view do not correspond to the true fields-of-view of the ACS 
and Suprim-CAM mounted on the HST and on the SUBARU 
telescope. For computational efficiency we limit the fields-of 
-view in the HST simulated observations to 120 arcsec, which 
is wide enough to contain the Einstein rings of our clusters. For 



the SUBARU simulations, the fields-of-view are defined such 
to correspond to the same comoving scale on the lens plane, i.e. 
4 hr x Mpc, for clusters gl and g51. For cluster g72, we simu- 
late a wider field-of-view, corresponding to 5 /r'Mpc comov- 
ing in order to include a large substructure in the observation. 

2.3. X-ray simulations 

The X-ray images of our simulated clusters are produced with 
the code X-ray MAp Simulator (XMAS). The software is pre- 
sented in Gardini et al. (2004) and Rasia et al. (2008), where a 
full description of the simulation pipeline can be found. It gen- 
erates synthetic event files which have the same format as real 
observations. In this work, which is not intended to exploit any 
X-ray calibration issue, we assume a constant response over 
the detector. In particular, the response matrix files and ancil- 
liary response file are those of the aimpoint of ACIS-S3 CCD 
on board of Chandra telescope. For each cluster projection we 
produce an X-ray image. The X-ray images have a field of view 
of 16 arcmin on a side, which, for the considered cosmology, 
corresponds to 2.5/r'Mpc at z = 0.2335 and to 2.97/T 1 Mpc 
at z = 0.297. The Chandra field-of-view is overlaid to the sur- 
face density maps of the nine clusters in Fig. 1 (green box). 
The background is not included a priori since Rasia et al. 2006 
showed that it does not induce any systematic on the bias. The 
spectral model used to generate the photons considers the con- 
tributions from the different metal species present in the simu- 
lation: C, N, O, Mg, Si, and Fe. The exposure time is 500 ksec. 
The images for all the cluster projections analyzed in this pa- 
per are shown in Fig. 4. The color-bar on the bottom allows to 
convert the color levels into counts per pixel. 

3. Analyses 

In this Sect, we describe the methods used for analyzing the 
previously outlined simulations. Firstly, we consider the mass 
estimates based on strong and weak lensing separately. Then 
we also discuss a non-parametric method which combines both 
the lensing regimes. Finally, we consider two methods for de- 
riving the cluster masses from the X-ray simulated data. 

3.1. Strong lensing 

The strong lensing analysis is performed by using the pub- 
lic software Lenstool (Kneib et al. 1993). This is very well 
developed tool for strong lensing parametric reconstructions, 
which allows to fit the observed strong lensing features in a 
cluster field through the combination of several mass compo- 
nents, each of which can be characterized by a density profile 
and by a projected shape (ellipticity and orientation). The code 
uses a bayesian approach to find the best-fit lens model and to 
estimate the errors on the free parameters. We refer the reader 
to the paper by Julio et al. (2007) for a more detailed descrip- 
tion of the many options available in this software. 

Lenstool allows the user to choose among several avail- 
able density profiles to describe the lens components. In this 
work, we use their implementation of the NFW profile (Golse 
et al. 2002) to model the main cluster halo. The functional form 
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Fig. 3. Left panel: Color-composite image of a simulated galaxy cluster, obtained by combining three SUBARU exposures of 2500s each in the 
B, V,/ bands. The field-of-view corresponds to ~ 450" x 450". Right panel: an HST/ACS composite image of the the central 100" x 100" of 
the same cluster. The image has been produced by combining mock observations with the filters F475W, F555W and F775W. 



of this profile is given in Eq. 2. Moreover, we add several sub- 
components, representing the contribution from the most mas- 
sive galaxies in the cluster. As shown in some previous studies, 
it is important to include the cluster members in the model, be- 
cause they can affect the positions and the magnifications of the 
strong lensing features (Meneghetti et al. 2003a; Meneghetti 
et al. 2007). These are modeled using Pseudo-Isothermal- 
Elliptical-Mass-Distributions (PIEMD) described by the fol- 
lowing density profile: 



PpiemdM - 



Po 



(l+r 2 /r c 2 ore )(l+r 2 /r c 2 u t) 



(12) 



We choose the PIEMD model because this is widely used for 
modeling the lensing properties of cluster galaxies in obser- 
vations (see e.g. Limousin et al. 2007; Riemer-S0rensen et al. 
2009; Donnarumma et al. 2009, for some recent references). 
As shown in the previous equation this profile is parametrized 
by a central density po, which is linked to the central velocity 
dispersion <r , and by two characteristic radii, namely the core 
radius r core and the cut-off radius r cut . To incorporate the galaxy 
population into the global lens model we use the same approach 
used by Limousin et al. (2007), who scale the parameters as a 
function of the luminosity (or equivalently of the mass). The 
scaling relations are given by 

L \'l 2 



core \Z*/ 



CT = CT, 



(13) 



where L* is a reference luminosity and the quantities r* ore , r* ut , 
and <t* are the corresponding core and cut radii, and central 
velocity dispersion, respectively. In our analysis, we assume 
a constant mass-to-light ratio for all the cluster members. Such 
an assumption is generally adopted when studying real clusters. 

Observationally, the galaxies to be included in the model 
should be selected as those laying in the cluster red sequence 
and being brighter than a given apparent luminosity (e.g. 
Limousin et al. 2007). Of course what matters for lensing is 
not the luminosity but the mass, which is assumed to be traced 
by the light. Indeed, the minimal luminosity should be inter- 
preted as a minimal mass. Working with simulations, we iden- 
tify the cluster galaxies using the SUBFIND code (Springel et al. 
2001) and then apply a selection based directly on the stellar 
mass. SUBFIND decomposes the cluster halo into a set of dis- 
joint substructures and then identifies each of them as a locally 
overdense region in the density field of the background halo. 

In our reconstructions, we include those galaxies which 
have stellar mass M slars > 10 9 hr l M Q and which are contained 
in a region of 500 /i^'kpc around the cluster center. This is typ- 
ically more than three times the size of the Einstein rings of 
the clusters in our sample. The orientation and the ellipticity of 
each galaxy are measured from the distribution of the star parti- 
cles belonging to it. Following this procedure, we typically end 
up with catalogs of several tenth of cluster members. 

The Brightest-Central-Galaxy (BCG) is included in the lens 
model by optimizing its parameters individually, rather than 
scaling them with the luminosity/mass. Since the BCG forms 
in the simulations in a strong cooling region, we assume it 
might have significantly different properties compared to the 
other cluster members. Thus, we prefer to treat it individually. 
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Fig. 4. X-ray maps of the three projections of the clusters gl (top panels), g51 (middle panels) and g72 (bottom panels) used in this work. The 
size of each map is 16 arcmin, which corresponds to 2.5/r'Mpc for g51 and to 2.97 /r'Mpc for gl and g72. In each panel, we display the 
masked regions. The dense cold blobs are encircled in green, while we indicate to the excluded central region in white. The color-bar on the 
bottom allows to convert colors into counts. 



Analogously, we use individual optimization with some other 
cluster members which lay particularly close to some multi- 
ple image systems. Indeed, their influence on the local lensing 
properties of the cluster requires to be carefully modeled. 

The total number of free parameters in the model depends 
on the complexity of the lens. Usually we consider a cluster- 
scale mass component, a galaxy-scale component to describe 



the BCG, and other galaxy-scale terms to incorporate the rele- 
vant cluster members. 

We distribute the sources behind the clusters such to have 
~ 3 - 7 strong lensing systems available for the optimization. 
For this condition to be satisfied, we randomly distribute few 
sources in a shell surrounding the lens caustics, enhancing the 
chances that they are strongly lensed. Then, we visually check 
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Fig. 5. Left panel: HST-like view over the central 120" x 120" of cluster gl (projected along the y axis). The light emission from the cluster 
galaxies has been removed to avoid confusion with the background multiple images. The identified systems are shown in yellow. The images ( 
belonging to the source n are indicated with the labels n.i. The inset on the left shows a zoom over the very central region of the image, where 
the central images 1.2, 3.4, 4.2, 5.2, and 6.4 are located. We also show in green the location of the cluster members which have been included 
in the lens reconstruction. Right panel: the result of the reconstruction obtained by using Lenstool with the lensing constraints shown in the left 
panel. The critical lines for sources at redshift 1.7 are shown in red and the corresponding caustics in yellow. We also overplot the true critical 
lines of the cluster in green. Note that the shape of the critical lines is very well reproduced, especially where the lensing constraints are tighter. 



if the multiple images belonging to each source are detectable 
in the simulation and retain those systems which are useful for 
the strong lensing analysis. The optimization is done using the 
bayesian method implemented in LENSTOOL with optimization 
rate 5 A = 0.1. We assume the uncertainty in the lensed image 
positions to be 07 = 0.3". 

In Fig. 5, we illustrate the reconstruction of cluster gl - y. 
The system has a massive galaxy at the center, labeled Gl, and 
another massive galaxy located ~ 45" west of the BCG (G2). 
Thus, we fit the lensing observables using three lens compo- 
nents, namely the main halo, the BCG and a secondary PIEMD 
clump coincident with the other massive galaxy. In this par- 
ticular case, adding additional cluster members does not af- 
fect the reconstruction. Among the sources, which were dis- 
tributed along the caustics of the input cluster lens, seven of 
them produce multiple image sytems detectable in this deep ex- 
posure (7500s) in the F775W filter. More precisely, two sources 
(source 3 and 6) produce five images, while the other sources 
are imaged into triplets. These are displayed in the left panel 
of Fig. 5. The bright knots in the multiple images of the same 
source are marked with points and identified with labels. The 
first digit corresponds to the source number, while the last in- 
dicates the multiplicity of the image to which the knots belong. 
For the sake of simplicity, the simulation is shown without in- 
cluding the light emission from the cluster members neither 
from other background sources that are not strongly lensed. 
The central image of source 2 cannot be detected by eye, thus 
it is not used in the reconstruction. The other central images of 



sources 1, 3, 4, 5, 6 lay behind the BCG but are detectable in 
the BCG subtracted frame. Using these observables, the recon- 
struction converges, finding a good fit to the lensing features 
ix 2 = 18 for 21 degrees of freedom). The best fit model consists 
of an NFW halo with concentration c = 10.57^^ ^ and scale ra- 
dius r s - 29.07^j° 8 ^ 8 arcsec (corresponding to 90^g 4 /i~'kpc). 
The galaxies Gl and G2 have velocity dispersions of 340+ jq 
km/s and 269^ km/s, respectively. To illustrate the result of 
the modeling, we show in the right panel of Fig. 5 the true and 
the reconstructed critical lines of the lens, assuming a source 
redshift of z s = 1.72 (source 3). These are displayed in green 
and red, respectively. Both the tangential and the radial criti- 
cal lines of the cluster are generally quite well reproduced by 
the model. The largest differences are on those portions of the 
critical lines along which fewer lensing constraints are present. 

The inner projected mass profile derived from the model, 
M(< R), is shown in Fig. 6 (dashed line), where we also show 
the true profile of the cluster acting as lens in this simulation 
(solid line). Since the model reproduces well the lens tangential 
critical line, it is not surprising that the model is very reliable at 
estimating the mass enclosed in the strong lensing region. The 
shaded area in the Fig. indicates the radial range of the multi- 
ple images, excluding the central images, which are located at 
R < 10 /r'kpc. The reconstruction reproduces well the true 
mass profile up to ~ 150 /r'kpc from the center, where the de- 
viation from the true mass profile is < 10%. At larger radii, the 
differences become significant. Thus, extrapolating the strong 
lensing model to distances where no strong lensing features are 
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Fig. 6. Results of the strong lensing analysis. The total projected mass 
profile of the inner region of cluster gl -y as recovered from the strong 
lensing mass reconstruction using LENST00L. The dashed line shows 
the result obtained by using seven multiple-image systems. The red 
three-dot-dashed line shows the mass profile if only three multiple- 
image systems are used. The blue dot-dashed line indicates the mass 
profile recovered by fitting all the 7 multiple-image systems, but as- 
suming that the all the central images are not detectable. The lensing 
constraints are shown in the left panel of Fig. 5. Finally, the true mass 
profile, as drawn from the particle distribution in the input cluster, is 
given by the solid line. The shaded region shows the radial range of 
the tangential strong lensing constraints. The bottom panel shows the 
ratios between the recovered mass profiles and the true mass profile. 

observed may result in largely incorrect mass estimates. This 
issue will be discussed in more detail in Sect. 4. 

In order to evaluate how the reliability of the model de- 
grades by reducing the number of constraints, we perform an- 
other reconstruction using only one system with five images 
(arising from the source 6) and two triplets (source 5 and 7). 
The final reconstruction does not differ significantly from the 
previous one. The projected mass profile for this new lens 
model is given by the red three-dot-dashed line in Fig. 6. This 
result shows that reliable reconstructions can be achieved even 
with a limited number of lensing constraints, if they are op- 
timally distributed across the cluster. We also attempted a re- 
construction by neglecting the presence of the central images 
(and using all the seven lensed systems). This is likely to be 
a realistic situation, since the central images are generally de- 
magnified and hidden behind the BCG, and thus difficult to de- 
tect. In this case, the mass enclosed by the strong lensing re- 
gion is again correctly estimated, but the reconstructed profile 
deviates more from the true one at small radii, as shown by the 
blue dot-dashed line in Fig. 6. In the following, the SL mod- 
els are constructed using also the central images, when they are 
detectable in the galaxy subtracted frames. 

3.2. Weak lensing 

The weak lensing measurements are done using the standard 
KSB method, proposed by Kaiser et al. (1995) and subse- 
quently extended by Luppino & Kaiser (1997) and by Hoekstra 
et al. (1998). Such method is now internally implemented 



in Skylens. The galaxy ellipticities are measured from the 
quadrupole moments of their surface brightness distributions, 
corrected for the PSF, and used to estimate the reduced shear 
under the assumption that the expectation value of the intrinsic 
source ellipticity vanishes (see Eq. A. 14). 

Selecting the galaxies with S /N > 10, we end up with cata- 
logs of galaxy ellipticities containing ~ 30 sources/sq. arcmin. 
The median redshift of these sources is z Sjt i Ue ~ 105. In the fol- 
lowing analysis we assume that all sources have the same red- 
shift of z s = 1. Furthermore, we assume that we can separate 
perfectly the population of background galaxies from the fore- 
ground cluster members. This is intentionally very optimistic, 
since we aim here at verifying the capabilities of several lensing 
methods to retrieve the cluster mass in the best possible condi- 
tions. The misidentification of cluster members as background 
galaxies leads to a dilution of the lensing signal, which causes 
to erroneous mass estimates (see e.g. Medezinski et al. 2007). 
We will address in more detail this issue in a forthcoming pa- 
per. Increasing the distance from the cluster center, the proba- 
bility that nearby substructures or additional mass clumps af- 
fect the mass estimates becomes higher. In this work, we do not 
include the effects of uncorrelated large-scale-structures (LSS) 
on the weak-lensing signal. The effects of the LSS on the weak 
lensing mass estimates have been discussed in detail in sev- 
eral other papers (Hoekstra 2001, 2003; White & Vale 2004). 
Uncorrelated LSS introduce a noise in the mass estimates, but 
the importance of matter along the line sight is fairly small for 
rich clusters at intermediate redshifts, like those in our sample, 
provided that the bulk of the sources are at high redshift com- 
pared to the cluster. Given that we are taking into account all 
the mass in cylinders of height 20 ft -1 Mpc in the lensing simu- 
lations, the effects of the correlated large-scale structure is par- 
tially included. Clowe et al. (2004) have studied the weak lens- 
ing signal of our same clusters (but in a pure dark matter ver- 
sion) in their cosmological environment. They find that includ- 
ing all the matter in a cylinder of height ~ 100/T 1 Mpc around 
the clusters only results in small scatter being added to the mea- 
sured cluster masses, compared to simulations where a cylinder 
of only ~ 10/i _l Mpc was used. These results are in agreement 
with the findings of Reblinsky & Bartelmann (1999). 

The cluster masses are derived using the following ap- 
proaches: 

NFW fit of the tangential shear profile: Assuming that the 
cluster is well described by an NFW density profile, we use 
the corresponding formula for the reduced shear to fit the az- 
imutally averaged profile of the tangential component of the 
reduced shear. For the NFW profile, the formulas for the radial 
profiles of the shear and of the convergence can be found in 
Bartelmann (1996) and in Meneghetti et al. (2003b). The tan- 
gential component of the reduced shear is given by 

g + = -Re[£e- 2 ^] , (14) 

where the angle <p specifies the direction from the galaxy cen- 
troid towards the center of the cluster, which we identify with 



M. Meneghetti et al.: Weighing simulated galaxy clusters using lensing and X-ray 



11 



+ 



0.3 



0.2 



0.1 



0.0 



-0.1 



1 1 1 1 


1 1 1 1 

9. » 
9. A 

NFW fit 


! | 1 1 ( 1 t i 

' ' ' ' ' 


1 t i ; 



100 200 300 400 500 600 
R [orcsec] 

Fig. 7. Radial profiles of the tangential and of the cross components 
of the reduced shear measured from the center of cluster gl — y. The 
dotted line shows the best fit NFW model. 



the position of the BCG. The cross component of the reduced 
shear is given by 



= -Im[ ge - 2 '*] . 



(15) 



If the distortion is caused by lensing, this component of the 
shear should be zero. 

In Fig. 7, we show the radial profiles of both the compo- 
nents of the shear for the cluster gl - y, measured out to large 
radii (~ 3 /i~'Mpc from the cluster center). The tangential com- 
ponent is well fitted by an NFW profile with c — 4.82 ± 0.64 
and r s = 0.307 + 0.048 /i^'Mpc. As expected in absence of sys- 
tematics, the cross component of the shear is consistent with 
zero. 

Aperture Mass Densitometry: The aperture mass densitome- 
try (Fahlman et al. 1994; Clowe et al. 1998) uses the fact that 
the shear can be related to a density contrast. More precisely, it 
can be shown that the following relation holds: 



((Ri) = k(< Ri) - k(R 2 <R< R mdX ) 



r 

Jr, 



d In R(y t ) 



2Rl 



f 



d\nR{y,) , (16) 



where k(< R\) is the mean convergence within a circular aper- 
ture of radius R\ and k(R 2 < R < R m ax) is the mean conver- 
gence in an annulus from R 2 to R llmx . This relation shows that 
the mean surface density within a circle can be derived from 
the tangential shear profile up to a constant. This is straightfor- 
wardly converted into a mass estimate: 



M(<Ri)=7rR 2 1 K(<R 1 )X ct (zi,Zs) 



(17) 



If wide field observations are available, R 2 and R mdx can be 
chosen to be large, so that the surface density in the annulus is 
negligible. Otherwise, setting the annulus term to zero, £ allows 
to estimate only a lower limit to the projected mass. 

Although the ^-statistic would not require any parameter- 
ization of the lensing signal to convert the shear into a mass 
estimate, a complication arises from the fact we do not directly 



measure the shear but the reduced shear. To convert the ob- 
served signal into an estimate of the shear, it is usually neces- 
sary to make some assumption on the shape of the convergence 
profile. In our analysis, we follow the method of Hoekstra 
(2007), who uses the convergence from the best fit NFW model 
to the tangential shear profile. We also use this approach to es- 
timate the mean surface density in the annulus. For our mass 
estimates R max varies between ~ 640 and ~ 800 arcsec, de- 
pending on the cluster. We set R 2 = 0.9 x R m!ix . 

3.3. Strong and weak lensing 

Finally, we consider a completely non parametric, two- 
dimensional mass reconstruction. The method used here is not 
based only on weak lensing. Instead, it combines both strong 
and weak lensing constraints to provide a map of the lens- 
ing potential. The method is fully described in Merten et al. 
(2009) (see also Bradac et al. 2005, for another similar al- 
gorithm). To obtain the underlying lensing potential if/ of the 
galaxy cluster the reconstruction algorithm performs a com- 
bined ^-minimisation, which consists of a weak and a strong- 
lensing term, in combination with an additional regularisation 
term: 



Rty). 



(18) 



The algorithm is grid-based and expresses the derivatives of 
the lensing potential by finite-differencing schemes. Starting 
from a coarse initial grid, the resolution is steadily increased 
until the final reconstruction is reached. At each iteration the 
reconstruction is regularised on the former iteration by the reg- 
ularisation function R(iff). This procedure results in a smooth 
reconstruction of the lensing potential and prevents the recon- 
struction from following noise patterns in the data. 

The weak lensing constraints are obtained by averaging 
over a certain number of background-galaxy ellipticities in ev- 
ery reconstruction pixel. The error is given by the statistical 
scatter in each pixel. Afterwards the reduced shear of the clus- 
ter is fitted: 



(€)■ 



z ^yw ) c -u (€) 



Z(z)y(<A) \ 
1-Z(z>wj. 



■(19) 



The sum over i and j runs over all reconstruction pixels and the 
cosmological-weight function Z(z) describes the redshift distri- 
bution of the sources 



Z(z)= 7 ^-H(z-z l ), 



(20) 



where D M and Di™ are the angular- diameter distances between 
observer and infinity and between lens and infinity respectively. 
Note that we have to solve the full y 1 -function since the recon- 
struction pixels might become correlated on a high grid resolu- 
tion. 

The strong lensing constraints are based on the estimated 
positions of the critical curves of the galaxy cluster. They are 
given by the observed arc positions and/or multiple images 
bracketing the critical lines. Assuming a reconstruction pixel 
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to be part of a critical curve, by summing over all these pixels 
with index k, we find: 



(d-z(zM<A)) 2 -|Z(z)rWI 2 ) 



(21) 



The error cr s is determined by the pixelisation of the grid. 
A final high resolution step is added, which is able to resolve 
the positions of the critical curve estimators in greater detail. 
Since there is no reliable weak lensing information on this res- 
olution, this step is embedded in the foregoing reconstruction 
by regularising on the former result in the cluster centre: 



*highresW =XtW)- 



R 



lowres 



(22) 



The mass analysis is then done with a convergence map, which 
is calculated from the reconstructed lensing potential of the 
cluster using Eq.A.5. 

In Fig. 8, we show the mass profiles of cluster gl — y ob- 
tained with the weak-lensing methods described above. All the 
methods perform well and provide consistent results for this 
particular cluster. The mass profiles deviate from the true one 
by less than 10% at all radii. The vertical lines show the posi- 
tion of /?2500> ^500) an d ^200) i- e - the estimated radii enclosing 
over-densities of 2500, 500, and 200 times the critical density 
of the universe, as derived from the NFW model that best fits 
the tangential shear profile. In the case of the two-dimensional 
mass reconstruction method combining strong and weak lens- 
ing, we estimate the errors on the masses by bootstrapping 24 
galaxy catalogs and by repeating the reconstruction with each 
of them. This is computationally very demanding, thus we re- 
construct the lensing potential a coarse grid of 32 x 32 pixels 
covering the whole cluster field (1280" x 1280"). 

3.4. X-ray 

Soft band [0.7-2] keV X-ray images are created from the event 
files. Using them we identify the regions of dense cold blobs 
that we mask in any further analysis. The masked regions are 
overlaid to the X-ray images for each cluster projection in 
Fig. 4. These bright point-like spots are mostly correlated to 
the cores of previously merged substructures. Further, we ex- 
clude an inner region of 60 kpc h~ l for g51 and 70 kpc hr x 
for gl and g72 (white circles in Fig. 4) to exclude the cen- 
tral region affected by the overcooling problem. From the soft 
band image, we finally produce the surface brightness profiles 
using 30 annuli spanning from 23 to 400 arcsec. The surface 
brightness profile of the cluster gl — y is given by the diamonds 
in Fig. 9. In the same radial range, we extract the spectra of 
10-15 annuli logarithmically spaced. We subsequently analyze 
the spectra in XSPEC (Arnaud 1996) using a single tempera- 
ture MEKAL model (Mewe et al. 1985; Liedahl et al. 1995) 
to obtain the temperature and the normalization. The values of 
redshift and hydrogen column density are set accordingly to the 
input simulation. The temperature profile of gl — y is shown in 
Fig. 10. 

The X-ray gas and total masses are recontructed using two 
different methods, here labelled backward wd forward method. 
In brief, starting from the observed X-ray surface brightness 
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Fig. 8. Results of the weak lensing analysis. Radial 2D-mass profiles 
of cluster gl - y, as obtained from the three different methods used 
in this work, namely the NFW fit to the shear profile (diamonds), the 
aperture mass densitometry (triangles), and the two-dimensional mass 
reconstruction combining weak and strong lensing (dotted line). The 
solid line shows the true mass profile. The vertical lines indicate the 
positions of i?2500, R500 an d ^2500 as derived from the NFW fit. The 
bottom panel shows the ratios between the mass profiles recovered 
from the lensing analysis and the true mass profile. The dotted line 
refers again to the SL+WL method, while the red and blue solid lines 
indicate the results for the NFW fit and for the aperture mass, respec- 
tively. 
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Fig. 9. Radial surface brightness profile of cluster gl - y, as derived 
from the analysis of the x-ray observation shown in upper middle 
panel of Fig. 4 (diamonds). The dotted and the dashed lines show the 
fits to the data obtained with the forward and with the backward meth- 
ods respectively. 



and the radially resolved spectroscopic temperature measure- 
ments and under the assumption of spherical geometry and hy- 
drostatic equilibrium between the ICM and the underlying dark 
matter potential, the backward puts constraints on a paramet- 
ric functional form of the total mass by deprojecting the ob- 
served quantities, while the forward makes direct use of the 
3-D model of the gas density and temperature profiles estimat- 
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Fig. 10. Radial temperature profile of cluster gl - y, as obtained from 
the spectral analysis of the X-ray observation (data points). The dotted 
and the dashed lines show the fits to the data obtained with the forward 
and with the backward methods respectively. See text for more details. 



ing the model parameters by projecting the profiles and fitting 
them to the observed ones. 

We describe here in detail the two methods. 

Backward method. A functional form of the total mass has to 
be assumed. From this, and moving backward, the profiles of 
the observed quantities are recovered. In this work, the NFW 
model in Eq. 2 is used to describe the density profile of the 
clusters, so that 



M tot (< r) = An r] p s f(x) , 
where 

/( O-Ind - v) 



1 + x 



(23) 



(24) 



and x = r/r s . 

The two free parameters (r s , C200) are constrained by mini- 
mizing ax 2 statistic defined as 



X 



2 _ \ 1 (^dataj ^Vnodelj) 



(25) 



where 7d ata are the either deprojected or observed temperature 
measurements obtained in the spectral analysis; T mo i e \ are the 
either 3-D or projected values recovered from the inversion of 
the hydrostatic equilibrium equation for a given gas density and 
total mass profiles; e T is the error on the spectral measurements. 
The gas density profile, n gas , is estimated from the geometrical 
deprojection (Fabian et al. 1981, Rriss et al. 1983, McLaughlin 
1999, Buote 2000, Ettori et al. 2002) of the measured X-ray sur- 
face brightness and, in the present analysis, is used to project 
J3D accordingly to the recipe in Mazzotta et al. (2004) to re- 
cover r mot jei- The latter values are then compared to the results 
of the XSPEC analysis using the^ 2 statistics (Eq. 25). The val- 
ues of are obtained from 



„ n gas M tot (< f) d (n gas X r 3D ) 
-Gfinip - 2 = , 



where G is the gravitational constant, m p is the proton mass and 
//=0.59 is the mean molecular weight in a.m.u. as adopted in 
the present simulations. Further detail on variations and appli- 
cations of the technique here described are presented in Ettori 
et al. (2002) and Morandi et al. (2007). The best fit surface 
brightness and temperature profiles obtained with this method 
for the cluster gl - y are given by the dashed lines in Figs. 9 
and 10, respectively. 

Forward method. By this technique, the total mass profile is 
recovered through a direct (forward) application of the hydro- 
static equilibrium equation (see Eq. 26) once a parametric form 
of the gas density and temperature profiles are estimated. In the 
following analysis, we adopt the approach by Vikhlinin et al. 
(2006) that follows similar techniques presented in, e.g., Lewis 
et al. (2003) and Pratt & Arnaud (2003), and extends the num- 
ber of parameters to increase the modeling freedom. In particu- 
lar, the formula describing the gas density is a y6-model modi- 
fied to accomodate a power-law behaviour in the center and the 
observed steepening of the surface brightness in the outskirts 
(Vikhlinin et al. 1999, Neumann 2005, Ettori & Balestra 2009) 
with the addition of a second y6-model to better reproduce the 
core profile: 

(r) _ Ni 1 (2J) 

N 2 

+ [l + Wr c2 ) 2 ] % 

a,j6i,/?2,'" J ,e, y,r c i, and r C 2 are all free parameters in the 
fit. The temperature profile is modeled by a 5 parameter 
-To, r t , a, b, c- function 

(r/r t T a 



T = T - 



Cool- 



(28) 



(26) 



For these simulations, where the cool cores are well inside 
the excluded inner regions (see Fig. 4), we consider = and 
Tcooi = L These profiles are then projected along the line of 
sight considering Mazzotta et al. (2004) prescription and the 
best-fit parameters determined from a x 1 minimization tech- 
nique by comparing the projected quantities with the observed 
ones. The best fit surface brightness and temperature profiles 
for cluster gl —y are given by the dotted lines in Figs. 9 and 10. 

To summarize, by using the observed X-ray surface bright- 
ness and temperature profiles these two methods provide: 
(backward method) starting from a given gas density profile, 
the gas mass M gas and the best-fit parameters (r s ,C2oo) from 
which a total mass profile and a 3D temperature profiles are re- 
covered; (forward method) for some given parametric forms of 
the gas and temperature profiles, their best-fit parameters from 
which M eas and M tot are estimated. In the following analysis, 
the results obtained from these two methods are compared to 
assess systematics in the X-ray analysis of the matter distribu- 
tion. 

In Fig. 1 1 we show the 3D-mass profiles of cluster gl — y, 
obtained with both the forward (triangles) and the backward 
(diamonds) methods. In the region between #2500 an d ^200 the 
two methods are consistent with each other and they under- 
estimate the true mass profile, given by the solid line, by 
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Fig. 11. Results of the X-ray analysis. Radial 3D-mass profiles of clus- 
ter gl - y, as obtained from the forward (triangles) and from the back- 
ward methods (diamonds) used in this work. The solid line shows the 
true mass profile. The vertical lines indicate the positions of ^2500, 
Rsoo and i?2500 as derived from the lensing (see Fig. 8, dashed lines) 
and from the X-ray backward analyses (dotted lines). 



Fig. 12. Ratio of the estimated to the true characteristic radii calcu- 
lated for each cluster projection. Squares, Stars, and Crosses refer to 
size derived using the weak lensing, the X-ray forward, and the X-ray 
backward method, respectively. The continuous line indicate where 
the true and the estimated size are equal. The Top, Middle, and Lower 
panel refer to R2500, and R200, respectively. 



~ 10 - 20%. Similar results are found for the other clusters 
in the sample, and will be discussed in detail in Sect. 4.3. The 
vertical lines in the figure indicate the estimated sizes of /?2500> 
R500, and /?2oo- In particular, the dashed and the dotted lines 
refer to the estimates based on the weak-lensing and on the X- 
ray forward analyses, respectively. Note that, due to the mass 
under-estimate in the X-ray case, the radii measured through 
lensing are typically larger than those measure from the X-ray 
methods. This is clearly shown in Fig. 12 where we report the 
ratio of the extimated to the true characteristic radii calculated 
for each cluster projection using the weak lensing signal or the 
X-ray analysis. It is worth noting even if the X-ray method al- 
ways underestimate the true radii, this bias seem to have a quite 
small scatter. On the contrary while the weak lensing is a less 
biased estimator, we find for it a much larger scatter. Just for 
conveniece, in the following analysis we compare the X-ray 
and the lensing masses at the same physical radii, which we 
choose to be /?2500, ^500, ^200 as derived from lensing. 

4. Results 

In this Sect, we discuss the results of the analyses outlined in 
the previous sections. We start with a discussion of the two- 
dimensional mass estimates obtained with the lensing tech- 
niques. Then, we consider the deprojection of the lensing pro- 
files and the X-ray 3D-mass estimates. Finally, we compare 
lensing and X-ray mass profiles and discuss how our results 
match the observations. 

4.1. Lensing 2D mass profiles 
4.1.1. Strong lensing masses 

In Fig. 13, we compare the true and the estimated strong- 
lensing masses of all clusters in our sample. The two- 
dimensional masses are measured at the limits of the strong 
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Fig. 13. The projected masses estimated through the strong lensing 
analysis vs the corresponding true masses of the lenses. The dotted 
lines correspond to M 2 d,sl = M tmt and to M 2 d,sl = M, ml , ± 10%. The 
masses are measured within a circle centered on the BCG and having 
a radius equal to the mean distance of the lensing constraints from the 
cluster center. 



lensing regions, given by the mean distance of the tangential 
images from the cluster center. The agreement between the es- 
timated and true masses is remarkable, showing that, in most 
cases, strong lensing methods based on parametric modeling 
are accurate at the level of few percent at predicting the pro- 
jected inner mass. The worst results are obtained for clusters 
g5 1 - z and g!2 - z, for which the offsets between the true and 
the estimated mass are ~ 15% and ~ 20%, respectively. These 
two clusters have complex morphologies, as shown in Fig. 1, 
being characterized by multiple mass components and double 
cores. For g51 - z the mass model combines 12 mass compo- 
nents, three of which are optimized individually. The cluster 
g72 - z is modeled with 78 components. In this projection the 
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Fig. 14. The relative difference between true and estimated projected 
masses as a function of the distance from the cluster centers. The re- 
sults are obtained by averaging over all the clusters in the sample. The 
diamonds and the triangles refer to the simulations including and ex- 
cluding the contribution of the BCG to the lensing signal (see text for 
more details). The errorbars show the scatter among all the reconstruc- 
tions. The regions probed by SL are typically smaller than ~ 30 - 40 
arcsec, thus the mass estimates at larger radii are extrapolations of the 
SL mass model to distances which are un-constrained by the data. 



cluster appears as a double cluster, with two main mass clumps 
separated by ~ 50". Each of these clumps is modeled with an 
NFW profile and has a massive star concentration at the center, 
which is modeled as a PIEMD. These four components are op- 
timized individually. Given the complexity of the systems and 
the large number of free parameters (24 for g72 - z), it is rea- 
sonable to expect less accurate mass estimates for these two 
clusters. 

Although the mass measurement is very precise at the posi- 
tion of the tangential critical lines, as anticipated in Sect. 3.1, 
the extrapolation of strong lensing models to larger radii such 
as 7^2500, ^5uo, or 7?2oo could lead to wrong mass estimates. 
In particular, we find that the results are very sensitive to 
the parameterization chosen to model the BCG, when we in- 
clude the central images in the optimization. Using the wrong 
model leads to biased mass estimates. In our simulations we 
adopted a PIEMD which is a broadly used profile for describ- 
ing the central galaxies in clusters. Such model, however, is 
not adequate to describe the profile of the BCGs in our simu- 
lations. While the PIEMD model is isothermal, i.e. the surface 
density profile decreases with the distance from the center as 
Spiemd(^) oc the true surface density distribution of the 
stars in our numerical simulations is steeper, i.e. Sbcg k R~ 1 ' 7 - 
This implies that, by fitting the lensing constraints as we do, we 
impose that the BCG is more spatially extended, and an addi- 
tional contribution to the central mass has to be provided by the 
dark matter halo. This causes to systematically over-estimate 
of the halo concentration and under-estimate of the scale ra- 
dius, as shown in Table 3. For this reason, the masses extrap- 
olated to large radii are systematically under-estimated. Such 
a dependence on the BCG mass has been reported recently by 
Donnarumma et al. (2009) modeling the cluster MS2137 (see 



also Comerford et al. 2006). The diamonds in Fig. 14 show the 
mean relative differences between true and estimated projected 
masses of the clusters at different distances from the centers. 
As the distance from the cluster center grows, the amplitude of 
the bias in the mass estimates increases, being of the order of 
~ 30% at a distance of two arcmins. The size of the errorbars 
reflects the scatter among the reconstructions. Note that also 
the scatter grows as a function of the distance from the center. 

In order to better illustrate how the results shown above de- 
pend on the stellar component of the lenses, we ran a new set 
of lensing simulations using only the dark matter halos of the 
clusters as deflectors. Then, we repeated the fitting procedure 
and we derived the mass profiles as done before, but without 
modeling the BCGs. In this case, even extrapolating the pro- 
files to radii much larger than the size of the strong lensing 
region, we find a much smaller disagreement between true and 
estimated masses. In this case, at large radii the strong lensing 
masses tend to be on average only slightly larger than the true 
masses (< 5%). This is shown by the triangles in Fig. 14. 

4.1 .2. Weak Lensing masses 

We discuss now the results obtained from the weak-lensing 
analysis of the clusters in our sample. 

The accuracy of the mass estimates depends on the mor- 
phology of the lenses and on their substructures. No big sub- 
structures are contained into the cylinders in the case of clus- 
ter gl. As shown in Fig. 1, the cluster appears regular in all 
the three projections. Instead, clusters g5l and g72 have more 
complex morphologies. In particular, as explained above, g72 
has a massive companion (~ 1.4 x 10 14 /z~'Mpc) located at 
~ 2.5 /z 'Mpc from the center. When projected along the z-axis 
of the simulation box the substructure is very close to the clus- 
ter center (~ 300/i~'kpc). For this reason, this mass clump has 
been included in the strong lensing model of g72 - z. Although 
it lays outside /?200 in the projections g!2 - x and g!2 - y, this 
substructure produces a significant shear which complicates the 
weak lensing mass measurements of the main cluster clump, as 
described below. 

In Fig. 15 we show the ratios of the estimated and of the true 
2D-masses of the clusters in the sample. The measurements 
have been made at /?2500 (diamonds), R^) (triangles), and /?200 
(squares). The different panels refer to the three methods im- 
plemented in this work, namely the NFW fit of the tangential 
shear profile (left panel), the aperture mass densitometry (cen- 
tral panel), and the SL+WL method (right panel). The reliabil- 
ity of the methods depends strongly on the morphology of the 
lenses. 

Two methods assume that the lensing signal is tangential 
to the center of the lens, which is assumed to coincide with 
the BCG. This implies that all shear is assumed to be pro- 
duced by a main cluster clump. These are the NFW fit and the 
aperture mass densitometry methods. For some lenses, like the 
three projections of gl, such approximation is rather well met. 
Indeed, this cluster is free of large substructures and the shear 
is dominated by the main cluster clump. For some other lenses, 
like the three projections of g72, this is not the case. In all three 
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Fig. 15. Comparison between the weak-lensing and the true 2D-masses of all the simulated clusters. From left to right, the panels refer to the 
methods based on the NFW fit of the shear profile, on the aperture mass densitometry, and on the non-parametric SL+WL reconstruction of the 
lensing potential, respectively. Shown are the ratios between the estimated and the true masses measured at three characteristic radii, namely 
^2500 (diamonds), R500 (triangles), and /?200 (squares), versus the cluster names. 



projections of this cluster the shear produced by the secondary 
mass clump contributes significantly to the total shear. As a 
result, in the projection along the x— and the y-axes, the tan- 
gential shear with respect to cluster center is smaller than it 
should be in absence of the secondary mass component. When 
fitting the shear profile with a single NFW model or converting 
the shear into a mass estimate using the aperture mass densit- 
ometry, this causes to under-estimate the mass. The amplitude 
of this effect grows with increasing distance from the cluster 
center, because the secondary mass clump is approached. In the 
case of g72 - x and g!2 -y, using the NFW fit method, the mass 
is underestimated by ~ 10 - 15% at fl 2 500 and by ~ 30 - 40% 
at /?20o- Similar results are found with the aperture mass den- 
sitometry. In the case of gll - z, where the secondary clump is 
located near the cluster center, the shear produced by the main 
cluster and by the substructure sum up. The resulting tangential 
shear mimics that of a lens with a larger mass and with an ex- 
tended core, because of the offset between the two mass clump 
centers. As a result, we find that the estimated M200 exceeds 
the true value by ~ 15 - 20%. 

In order to validate our interpretation of these results, we 
perform the following tests. First, we place two NFW halos 
of mass Mi = 6.8 x 10 14 h~ l M Q and M 2 = 1.4 x 10 14 at a 
distance of 2.5 hr l Mpc, and we calculate the shear field pro- 
duced by these two mass clumps. Then, we sample the shear 
field at random galaxy positions and add the noise due to the 
intrinsic galaxy ellipticities, following the method outlined in 
Maturi et al. (2005). Finally we measure the tangential shear 
profile with respect to the center of the most massive halo and 
we fit it using a single NFW profile, over a radial range similar 
to that used in our numerical simulations. We repeat the exper- 
iment after removing the least massive halo. We find that, if the 
shear produced by the second halo is included in the simula- 
tion, the masses at /?250Q> ^500> an d ^200 are under-estimated, 
with respect to the simulation including only the most massive 
halo, by 4%, 21%, and 30%, respectively. This is qualitatively 
in agreement with the results for g!2 - x and g!2 - y. Note 
that in both cases, the tangential shear profile is well fitted by 
an NFW profile, i.e. with reduced ;f 2 < 1 . A a second test, we 
place the least massive halo at a distance of 300 /i _1 kpc from 
the main halo. The resulting tangential shear profile is com- 
pared with that produced by a single mass component of total 



mass M = Mi + M 2 . By fitting with NFW models, we find that, 
in the two-halo case, the masses at /?2500> ^500> an d ^200 are 
over-estimated, by 3%, 17%, and 25%, respectively, with the 
respect to the single halo case. This is also in good agreement 
with our results for g72 - z. 

For the remaining clusters, the mass estimates obtained 
with these two methods are more accurate. Typically, the esti- 
mated masses deviate by < 15% from the true masses. Even in 
the case of g5 1, despite the presence of few substructures in the 
cluster surroundings, the mass estimates are in good agreement 
with the input masses. The reason is that the above mentioned 
substructure are less massive than in the case of g72 and the 
shear is dominated by the cluster halo. 

The non-parametric SL+WL method does not require 
any assumption on the symmetry of the lensing signal. 
Substructures are included in the mass model by construction. 
For this reason, the SL+WL method can recover the input mass 
with good precision even in the case of morphologically dis- 
turbed clusters. We find that the deviations between estimated 
and true masses are typically below the 10% level. 

Note that, for all the three methods the scatter between es- 
timated and true masses grows as a function of the radius, in 
agreement with Okabe et al. (2009). 

4.2. Lensing 3D masses 

As stated several times above, lensing measures the total mass 
contained in cylinders and projected on the sky. To convert the 
2D-mass into 3D-mass estimates, deprojection needs to be im- 
plemented. This requires to make some assumptions on the 
shape of the clusters and on their three-dimensional density 
profile. We assume here that clusters are spherical and that their 
density profile is well described by the NFW model. 

De-projection is done differently for the three methods in- 
vestigated here. For the NFW fit method, we use the NFW pa- 
rameters obtained from the fit of the tangential shear profile 
to calculate the 3D-mass profile of the lenses. For the aper- 
ture mass densitometry and for the SL+WL methods, we fit the 
2D-mass profiles with projected NFW models and we use the 
best fit parameters to derive the 3D-mass profiles. These 3D- 
profiles are shown in Fig. B.l for the individual cluster projec- 
tions. 
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Fig. 16. As in Fig. 15, but comparing the estimated and true 3D-masses. 



Similarly to Fig. 15, we show in Fig. 16 the ratios between 
estimated and true 3D-masses at three over-density radii. It is 
clear that in 3D the scatter between the estimates and the in- 
put masses is significantly larger than in 2D. This is due to two 
main reasons. First, lenses are triaxial, while we are assuming 
spherical symmetry during the de-projection. Depending on the 
degree of triaxiality and on the orientation of the clusters with 
respect to the line of sight, 3D-masses may result to be over- 
or under-estimated. In particular, we find that, in the cases of 
good alignment (i.e. small angles) between the major axis of 
the cluster and the projection axis, the lensing masses tend to 
be systematically larger than the true masses, while the oppo- 
site occurs in those cases where the major axis is nearly per- 
pendicular to the line of sight. This is shown in Fig. 17, where 
the lensing masses are derived with the SL+WL method, which 
even in 3D seems to provide the most accurate mass estimates. 
Given that the masses of ellipsoids and spheres with the same 
azimuthal density profile tend to converge at large distances 
from their centers, the effect is strongest for M2500 and for M500, 
and mildest for M2oo- However, even at R200, the analysis of our 
sample shows that the scatter due to triaxiality is of the order 
of ~ 20%. Similar results are found by Corless & King (2007). 

The second factor which makes the 3D lensing mass esti- 
mates so noisy is the presence of substructures along the line 
of sight. Since their distance from the lens plane is unknown, 
the 3D-mass estimates can be severely affected by these mass 
clumps, especially if they are located close to the cluster core 
in projection. The high ratio between the estimated and the true 
mass of g!2-z is in large part due to the presence of the massive 
sub-clump previously mentioned in the paper. This accounts for 
~ 15% of the total cluster mass, but its erroneous inclusion in 
the central 300/i~'kpc significantly affects the mass estimates 
at small radii. 

Apart from this particular cluster projection, in the cases of 
systems without large substructures along the line of sight, 3D- 
lensing masses are affected by an intrinsic uncertainty due to 
triaxiality, which we estimate to be of the order of 20% at /?2oo- 
Unfortunately, at smaller radii where the lensing measurements 
of the 2D masses would be more robust, the scatter becomes 
larger, being of the order of 50%. 

Recently Okabe et al. (2009), by studying the weak lensing 
signal of 30 galaxy clusters observed with the SUBARU tele- 
scope, found that the mean ratios between 3D- and 2D-masses 
at R vn and R 50 o are 1.34 + 0.17 and 1.40 + 0.10, respectively, 



1.4 



I 1.2 



1.0 



0.8 



g72-z 
g72-z 

g72-z 

□ 9 51 " 
A 951- 
O 951- 



M 2500 


O 


M 500 


A 


M 200 


□ 



gl-y 



g1- 

9'- 



□ 

A 
O 

g51 



g51- 
g51- 
g51- 



20 40 60 80 100 
Angle with the proj. axis [deg] 



120 



Fig. 17. Ratio between estimated and true lensing masses as a function 
of the angle between the major axis of the cluster inertia ellipsoid and 
the axis along which the mass distribution is projected. The results 
are shown for the lensing masses obtained with the SL+WL method. 
Squares, triangles and diamonds indicate the mass measurements at 
#200, #500, and R 25 oo, respectively. 



where is the virial radius. They derive the 2D-masses using 
the aperture mass densitometry method, while the 3D-masses 
are obtained from the NFW fits of the shear profiles. Their sam- 
ple spans a range of masses which is much wider than that 
covered by our sample. From their Fig. 8, we can estimate 
that, limiting the analysis to masses M vir > 8 x 10 14 hr l M Q 
and M500 > 4 x 10 14 hr 1 M B , the ratios are smaller (~ 
1.06 and ~ 1.37 at R v i r and at R500, respectively). Averaging 
over our sample, we find M\ 



' vir,2D 



/ M S> = L14 ± 009 and 



MZ,,JM?™„ = 1.29 + 0.08 



500.2D' 500.3D 

These are quite in a good agreement with the the ratios be- 
tween the true 2D- and 3D-masses, which are M,'™ n /M^™ e 3D = 



1.13 and M^ 2D /M^ 3D 



vir,2D' 

1.34. From these results, we can 



deduce that a) on average, the ratios between 2D- and 3D- 
masses are well recovered, despite triaxiality and substructures 
affecting individual mass estimates, and b) the agreement be- 
tween simulations and observations is an indication that the 
density profiles of real and simulated clusters are, within the 
errors, compatible with each other. 
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Fig. 18. Top panels: comparison between the X-ray and the true 3D-masses of all the simulated clusters. The left and the right panels refer 
to the forward and to the backward methods, respectively. Shown are the ratios between the estimated and the true masses measured at three 
characteristic radii, namely ^2500 (diamonds), R500 (triangles), and W200 (squares), as a function of the cluster name. Bottom panels: as for the 
upper panels, but comparing the gas masses derived from the X-ray analysis to the true gas masses of the clusters. 



4.3. X-ray masses 

In Fig 18 we compare the true and estimated total and gas 
masses for all the cluster in our sample. The left and the right 
panels refer to the forward and to the backward methods de- 
scribed in Sect. 3.4. The top panels show the ratios between 
the estimated and the true masses measured at three character- 
istic radii, namely /?250o (diamonds), Rsqq (triangles), and 7?2oo 
(squares), as a function of the cluster name. A comparison of 
the mass profiles at any other radius is shown in Fig. B.l. The 
forward and backward methods give results which, at the three 
different characteristic radii, are consistent within 2<x for all the 
clusters. The reconstruction of the X-ray masses through the 
two methods can differ because, while the backward method 
assumes a functional form of the total density profile impos- 
ing the gas temperature profile for given gas density profile, 
the forward method reproduces the gas density and temparture 
with very flexible models that, in particular in the outer regions, 
can induce not stable extrapolations. For the forward method, 
the confidence intervals on mass profiles have been estimated 
by minimising the distance between the projected models and a 
set of 100 random realisations of observed profiles. These pro- 



files have been obtained assuming Gaussian statistics around 
observed values of the surface brightness and temperature pro- 
files with the constraint of rejecting realisations leading to non- 
physical solution of non-monotonically increasing mass pro- 
files. For the backward method, errors on M gas x are estimated 
by propagation of the errors on the gas density profile as ob- 
tained from the direct deprojection of the surface brightness, 
reconstructed 100 times by considering the Gaussian error in 
each radial bin. The error on Afjo.x is obtained from the prop- 
agation of the errors on the concentration parameter and scale 
radius estimated through Eq. 25. Similarly to what shown in 
other work (see e.g. Rasia et al. 2006, 2008; Nagai et al. 2007, 
and references therein) we find that, overall the X-ray mass es- 
timates are biased low by 5%-20%. More quantitatively we find 
that, on average, both methods underestimate the true mass by 
10% with a scatter of 6% almost independently of the char- 
acteristic radius considered. The most deviant mass estimate 
comes only from the z projection of cluster g72. This is due 
again to the fact that this cluster is composed by two structures 
that line-up along the line of sight (see Fig 1). When we fit 
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the surface brightness profile this appears to be steeper in the 
center and shallower at the characteristic radius. 

In the bottom panels of Fig. 18 we compare the gas masses 
derived from the X-ray analysis to the true gas masses of the 
clusters. Again, we find that the two methods give results con- 
sistent within 2cr errors. Furthermore, in agreement with other 
work, we find that the total gas mass is well reconstructed 
(Nagai et al. 2007). More quantitatively we find that, on av- 
erage, the gas mass is recovered to better than 1% with a 3% 
scatter at /? 2 5oo and 7% with a 3% scatter at R 500 and R 2 qq. It is 
important to say that the trend to slightly overestimate the gas 
mass at large radii is due to the fact that both X-ray methods 
tend to slightly overestimate the cluster density profiles in the 
outskirts due to a non perfect azimuthal symmetry of the clus- 
ter surface brightness at larger radii. The most deviant gas mass 
estimates is given by the z projection of cluster g72. Again this 
is due to the same effect explained above for the mass. 

4.4. NFW fit parameters 

The best-fit NFW parameters for all the clusters in the sam- 
ple, obtained both through the lensing and X-ray analyses in 
the radial ranges listed in Table 2, are summarized in Table 3. 
In column 4 of the table we report the best fit concentrations 
and scale radii obtained for the main cluster halo using the 
strong-lensing constraints. These results are obtained by in- 
cluding also the central images in the strong lensing modeling. 
As discussed in Sect. 4.1.1, this leads to systematically over- 
estimate the concentration and underestimate the scale radius, 
being the PIEMD model used in the SL modeling inadequate 
to describe the distribution of the stars in the simulations. In 
particular, we note that, in some cases, the concentrations are 
off by a factor of ~ 3 with respect to the true values reported 
in the second column. On average, the concentrations derived 
from SL alone are almost 85% larger than the true concentra- 
tions of the DM-only profiles. Conversely, the scale radii are 
almost ~ 50% smaller. We remind that the SL fits refer to the 
cluster halos, while the stars are modeled apart.. Such bias is 
not present in the simulations without stars, where the aver- 
age ratio between estimated and true concentrations is around 
unity. 

The weak lensing best fit parameters, obtained from the fit 
of the shear profiles, are given in column 5. The concentrations 
tend to be slightly larger than those of the dark-matter distri- 
butions in the input models ((c WL /c™ only ) = 1.17) and the 
scale radii are on average ~ 10% smaller that the input values 
((r* L /r™° nly ) = 0.93). Similar results are found combining 
SL and WL (columns 6). In this case the average ratio between 
estimated and true concentrations is (c SL+WL /c™ only ) = 1.15 
(for the scale radius we find <rf L+WL /r™° nly > = 0.93). The 
tendency to over-estimate more the concentration when adding 
the SL constraints is caused by the large contribution of the 
stellar and gas masses within the inner 100 /r'kpc, as shown 
in Fig. 2. We remind that in the cases of the WL and WL+SL 
methods the fits are done over the total projected mass pro- 
files, thus without distinguishing the dark-matter component 
from the stellar and gas masses. To support this interpretation, 



we note that the largest discrepancies with the true concentra- 
tions arise for the cluster g5l, which is characterized by an ex- 
tended strong over-cooling region in the center, where the to- 
tal density profile steepens compared to the DM only profile. 
This mimics a larger concentration, being the mass profile still 
compatible with an NFW model. Fitting the total density pro- 
files of the inputs clusters indeed leads to higher concentrations 
and smaller scale radii, as reported in column 3, which are in 
a much better agreement with the results of the lensing fits. As 
highlighted in the previous sections, the WL fits of some sys- 
tems, like the three projections of g72, can be strongly biased 
assuming a single mass component. Indeed, their substructures 
need to be properly modeled when deriving the cluster mass 
from the shear signal. The SL+WL method provides a better 
chance to measure the density profile, at least in some cases, 
like g72 - x and g!2 - y, where the large substructure laying at 
the edge of the cluster virial region is by construction properly 
modeled. For g!2 - z, where the substructure is near the cluster 
core, a correct fit of the projected mass distribution would re- 
quire multiple mass components also for the SL+WL method. 

Comparable estimates of c and r s are obtained with the 
two X-ray techniques. Even the X-ray analysis probes the to- 
tal mass distribution, including the gas and the stars. We note 
that, using the Backward method, the concentration and the 
scale radius are obtained as best-fit parameters that minimize 
the reconstructed temperature profile with respect to the ob- 
served values, whereas using the Forward method we have es- 
timated them with a-posteriori fit with a NFW functional form 
performed in the radial range 0.7 - 1.4 /T 1 Mpc on the mass 
model obtained by appling the hydrostatic equilibrium equa- 
tion on the gas density and temperature models described in 
Eqs. (27) and (28), respectively. Overall, the deviations from 
the true estimates go in the same direction: in gl, c is underes- 
timated consistently by a factor 0.7 - 0.8, with a corresponding 
overestimate of r s up to a factor 1.45, if compared with the true 
concentrations obtained by fitting the input total mass distri- 
butions of the clusters; the same considerations apply to g72, 
apart from g72 - z, where the concentration (scale radius) is 
estimated higher (lower) than c l °^f (ff^) by about 15 per cent 
as consequence of the alignment along the line of sight of the 
two main clumps; more critical is the case of §51, where both 
the X-ray methods provide a measure of the concentration that 
is twice c^ 1 , due again to the large contribution of the cool 
substructures in the central regions, as also shown from the 
overestimate of the total mass obtained through X-ray analy- 
sis within /?2500 ( see Fig. 18). Although we attempted to mask 
it in the X-ray analysis, the over-cooling region in this cluster is 
very extended, thus it is still affecting the mass reconstruction 
(see Fig. B.l in Appendix B). On average the X-ray concen- 
trations are ~ 10% and ~ 20% smaller than c™ only and c™f, 
respectively, if the three projections of g5l are not included. 



4.5. Lensing vs. X-ray 3D mass profiles 

Finally, we attempt a comparison between the lensing and the 
X-ray mass estimates. 
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Table 2. Radial ranges used for fitting NFW profiles to the lensing and X-ray data. All radii are expressed in units of r 2 m reported in Table 1. 



cluster min-max SL min-max WL min-max SL+WL min-max X-ray 



gl- 


X 


0.002-0.08 


0.16-1.00 


0.03-1.00 


0.45-0.90 


gl- 


y 


0.002-0.09 


0.16-1.00 


0.03-1.00 


0.45-0.90 


gl- 


z 


0.002-0.06 


0.16-1.00 


0.03-1.00 


0.45-0.90 


g51 


- X 


0.003-0.07 


0.20-1.00 


0.05-1.00 


0.50-1.00 


g51 


-y 


0.002-0.05 


0.20-1.00 


0.05-1.00 


0.50-1.00 


g51 


- z 


0.002-0.12 


0.20-1.00 


0.05-1.00 


0.50-1.00 


g72 


- X 


0.006-0.05 


0.20-1.00 


0.05-1.00 


0.54-1.08 


gV2 


-y 


0.006-0.05 


0.20-1.00 


0.05-1.00 


0.54-1.08 


g72 


- z 


0.002-0.05 


0.20-1.00 


0.05-1.00 


0.54-1.08 



Table 3. The NFW concentrations and scale radii (upper and lower part of the Table, respectively) resulting from the strong-lensing (column 
4), weak-lensing (column 5), strong+weak lensing (column 6), and X-ray analyses (column 7-8 for the forward and for the backward methods) 
of the clusters in our sample. The weak-lensing estimates are obtained by fitting the shear profile with an NFW model. In columns 2 and 3, 
we quote the true concentrations obtained by fitting the DM-only and the total density profiles of the three clusters in the radial range between 
10 h' 1 kpc and r2oo- Below each column, we report the mean ratios of the estimated and the true parameters. These are calculated both including 
and neglecting the three projections of g51, which strongly bias the mean concentrations and scale radii derived from the X-ray analyses. The 
numbers in parentheses are the corresponding r.m.s. values. 



cluster 


DMonly 
C true 


total 
L tme 


C SL 


C WL 


C SL+WL 


^X.forw 


^.X.back 



gl-x 4.62 5.38 H.31° 19 59 6.21 ± 0.87 6.59 ± 0.32 3.58° 006 -^_ 028 

gl-y 4.62 5.38 10.57 2 f 2 81 4.82 ± 0.64 5.39 ±0.23 4.51° 08 08 4.05 - 40 25 

gl-z 4.62 5.38 6.92 2 f 12 4.44 ± 0.62 4.54 ± 0.23 4.06° 07 07 3.62 022 , 

g51 -x 5.37 7.20 3.36°f n 5.83 ± 0.96 6.10 ±0.29 U.15°f 20 10.76°f 65 

g51-y 5.37 7.20 9.72° 6 ' 8 61 5.13 ± 1.05 6.35 ± 0.42 12.39 a22 2 10.80°f 56 

g51-z 5.37 7.20 8.78°^ 6.06 ± 0.93 7.41 ±0.32 H-04° 2 ' 22 H-30 f 62 

g72-x 3.99 4.22 6.46°f 54 4.17 ± 1.54 3.88 ± 0.30 3.31° 05 05 3.22°° 6 02 

g72-y 3.99 4.22 7.26 1 /f 14 7.91 ± 2.86 4.17 ±0.31 3.48° 05 05 3.29°° 6 02 

g72-z 3.99 4.22 11.3 9°f 85 4.19 ± 0.62 4.50 ± 0.24 4.92 08 08 4.49 03 ' 8 29 



1 1 DMonly 
C l C true 
/total 
W L true 




1.84 (0.60) 
1.59 (0.59) 


1.17 (0.30) 
1.01 (0.33) 


1.15 (0.15) 
0.98 (0.11) 


1.33 (0.59) 
1.09 (0.38) 


1.25 (0.56) 
1.03 (0.36) 


/ DMonly 
c / c tn.e 
r /^total 
W L tnie 


no g51 
no g51 


2.09 (0.48) 
1.88 (0.45) 


1.23 (0.35) 
1.12(0.35) 


1.12(0.15) 
1.00 (0.12) 


0.92 (0.15) 
0.83 (0.15) 


0.86 (0.12) 
0.78 (0.13) 




DMonly 

r 

s,true 


..total „SL 
' s,true ' s 


„WL 

s 


SL+WL 

s 


X,forw 
' s 


^X.back 



gl-x 0.310 0.278 0.097° 00 7 02 0.229 ± 0.037 0.221 ±0.012 0.408° 00 7 7 0.409 038 7 

gl-y 0.310 0.278 0.090° °' 6 4 06 0.307 ± 0.048 0.276 ±0.013 0.315° 00 6 06 0.360°f 6 3f> 

gl-z 0.310 0.278 0.152%% 0.317 ± 0.053 0.320 ± 0.017 0.351%™, 0.403° 04 3 27 

g51 -x 0.241 0.189 0.385^3 0.242 ± 0.045 0.235 ± 0.012 0.120° 00 2 2 0.124° °° 9 06 

g51-y 0.241 0.189 0.099 00 8 2 0.246 ± 0.057 0.206 ±0.014 0.107° 00 2 2 0.124 °° 8 7 

g51-z 0.241 0.189 0.141 00 2 3 0.244 ± 0.043 0.204 ±0.010 0.120 00 3 03 0.116 ™ 7 

g72-x 0.299 0.299 0.087° ' 32 9 0.262 ± 0.105 0.306 ± 0.025 0.366° 00 6 06 0.386° *' 07 

g72-y 0.299 0.299 0.094° 029 27 0.135 ± 0.053 0.287 ± 0.023 0.355 00 5 06 0.386 00 2 8 

g72-z 0.299 0.299 0.059 ' 08 0.367 ± 0.062 0.324 ±0.018 0.223° 00 4 04 0.251° 01 9 22 



, DMonly 
r s/ r s,true 
„ /rtotal 
's/' s ,true 


0.49 (0.40) 
0.59 (0.53) 


0.93 (0.20) 
1.05 (0.27) 


0.93 (0.11) 
1.04 (0.12) 


0.89 (0.33) 
0.98 (0.32) 


0.97 (0.36) 
1.06 (0.34) 


, DMonly c , 
''s/'- s ,true n °g 51 


0.31 (0.09) 
0.34 (0.11) 


0.88 (0.24) 
0.94 (0.26) 


0.95 (0.12) 
1.001 (0.11) 


1.10(0.18) 
1.17 (0.22) 


1.20(0.17) 
1.27(0.21) 



First of all, given the results discussed above, it is not 
surprising that lensing and X-ray mass estimates in individ- 
ual galaxy clusters can differ by up to 100%. This is clear in 
Fig. 19, where the lensing masses are shown as a function of 
their X-ray equivalents. The lensing masses are obtained with 
the SL+WL method. The X-ray masses are obtained with the 
forward method. The discrepancies in the simulated sample are 
consistent with those observed in several galaxy clusters. As an 
example, we over-plot with asterisks the mass measurements 



of a sample of 18 galaxy clusters reported by Mahdavi et al. 
(2008) (M08 in the following discussion). 

Several recent studies seem to agree on the fact that the 
lensing masses are on average larger than the X-ray masses. 
For example, analyzing a sample of 19 clusters observed with 
SUBARU and XMM-Newton, Zhang et al. (2008) find that 
-^5oo,wl/'^500,x = 1.09 ± 0.08. M08 report similar results. In 
particular, they find a trend in the ratio between lensing and 
X-ray masses as a function of the overdensity radius. While 
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X-ray and the lensing masses are consistent with each other 
at #2500 (M 2 5oo,x/M 2 5oo,wl = 1.03 ± 0.07), their mean ratio 
becomes 0.78 + 0.09 at /?5oo- Correcting for excess correlated 
structure outside the virial radius, they find M2soo,x/M25oo,wx = 
1.06 + 0.07 and M 5 oo,x/M 5 oo,wl = 0.85 ± 0.10. Very re- 
cently, Zhang et al. (2009) (Z09 in the following), studying 
a sample of 12 clusters divided into relaxed and unrelaxed 
on the basis of their X-ray morphology, also find that the 
weak lensing masses exceed the X-ray masses at large radii 
(or low over-densities), while they are comparable at small 
radii (or high over-densities). However, they measure a shal- 
lower but still significant radial evolution of Mx/Ml com- 
pared to M08. For the sub-sample of relaxed clusters, they 
find M 2 50o,x/M 2 50o,wl = 1.04 ± 0.08 and M 5 oo,x/M 50 o,wl = 
0.91 ± 0.06 on average, while for the whole sample of objects 
they find M 25 oo,x/M 25 oo,wL = 0.97+0.07 and M 5 oo,x/M 50 o,wl = 
0.94 + 0.05. Such trend is interpreted as an indication of lack 
of hydrostatic equilibrium in galaxy clusters (Churazov et al. 
2008). 

We repeat here the analysis of M08 and Z09 using our sim- 
ulated clusters. Following the terminology of M08, we indicate 
with «a the ratio between the X-ray and 3D-lensing masses for 
a particular overdensity A. We measure a a for A = 2500, 1000, 
500, and 200. Consistently with M08 and Z09, we estimate a A 
by minimizing ax 2 statistic defined as 



X 



s 



(Af A , x - QaMax) 2 

°"L + fl A°AZ 



(29) 



where o"a,x and <Xax are the errors on the X-ray and lensing 
masses corresponding to the overdensity A. The errors on a& 
are estimated by by locating the values at which x 2 ~Xmia - 1. 
which correspond to the 68% condence interval. For this anal- 
ysis, we use the lensing masses obtained from the SL+WL 
method and the X-ray masses obtained from both the forward 
and the backward methods. Note that the lensing masses in 
M08 and Z09 are obtained with the NFW fit method. We opt 
for the SL+WL method because our sample is limited and a 
significant fraction of our clusters has significant substructures. 
Under these conditions the efficiency of the NFW fit method to 
obtain reliable mass estimates is limited, as we discussed in 
Sect. 4.1.2. The results are shown by the squares in Fig. 20, 
where we plot «a as a function of A. For comparison, we over- 
plot the data-points taken from Fig. 4 of M08. We also plot the 
results obtained by Z09 for their sub-sample of relaxed clus- 
ters.. We find that the ratios between X-ray and lensing masses 
are below unity at all over-density radii. For A = 2500, we find 
a a = 0.9*q q3 and a & ~ 0.91q^ using the masses derived from 
the forward and from the backward methods, respectively. For 
A = 500 such ratio becomes «a = 0.88^2 ™ using the forward 
method to measure the X-ray mass, and using the backward 
method it is «a = 0.87*2^. We do not detect a strong ra- 
dial trend as reported in M08 and Z09. On the contrary, the 
ratios between X-ray and lensing masses decline very gently as 
a function of the overdensity radius. . 

The important question which arises now is if the ratio be- 
tween lensing and X-ray masses is a robust indicator of the lack 
of hydrostatic equilibrium in the simulated clusters. This can 
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Fig. 19. Lensing vs X-ray masses within /?2soo (diamonds), ^ 50 o (tri- 
angles), and /?200 (squares). The asterisks show the mass estimates 
published by M08 for a sample of 18 galaxy clusters. The dashed lines 
correspond to M L = [0.5, 1.0, 1.5,2.0] x M x . 



be easily checked by measuring the mass profile of the simu- 
lated clusters using the hydrostatic equilibrium approximation. 
For doing this, we plug the gas density and the temperature 
profiles measured from the particle distributions into the hy- 
drostatic equilibrium equation and we derive the mass Mheq- 
This mass should be recovered by the X-ray analysis, i.e. Mx 
is a measurement of Mheq- This mass is compared to the true 
mass M true as in Rasia et al. (2004). By considering the three 
clusters, we obtain the median profile of the ratio between hy- 
drostatic equilibrium and total mass, shown by the dashed line 
in Fig. 20. For comparison, we indicate with asterisks the ra- 
tios between X-ray and true masses at each over-density radius. 
In these simulations, the median MHEQ/M tme is ~ 0.9 between 
/? 2 50o and /?iooo and then it decreases to 0.84 and 0.8 at R^qq and 
#200- The X-ray masses are compatible (within the 68% confi- 
dence limit) with Mheq- The mass derived with the backward 
method declines more steeply than that obtained with the for- 
ward method as A decreases, and seems to follow more closely 
the behaviour of M H eq, although the differences are not sig- 
nificant given the error-bars. Mx/Ml is also compatible with 
MHEQ/M true . This is true only at the 2<x significance level for 
A < 500. Thus, we conclude that lensing and X-ray obser- 
vations and their comparison to the simulations can provide 
important information on the physics of the gas in galaxy clus- 
ters. Note that, although we have considered a small sample 
of clusters, averaging over several objects significantly reduces 
the impact of triaxiality. Indeed, at all overdiensities, Mx/Ml 
is very similar to Mx/M tlue , indicating that Ml ~ M^e on av- 
erage. 

In Table 4 we report the mean, the r.m.s., the median, the 
first and the third quartiles, and the minimal and maximal val- 
ues of the distributions of several mass ratios. Not only we 
compare 3D vs 3D masses, but we also compare 3D vs 2D 
mass estimates and vice-versa. We note that triaxiality and sub- 
structure affect the 2D X-ray mass estimates consistently to the 
3D lensing mass estimates. For example the r.m.s. of the ratio 
M 2 d,x/M 2 d,l is of order ~ 17 — 18%, similar to the r.m.s. of 
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Fig. 20. Ratio between X-ray and lensing masses as a function of the 
overdensity A (squares). The results are shown for the X-ray masses 
obtained with the forward (red) and with the backward (blue) methods 
and for the lensing masses obtained with the SL+WL method. For 
comparison, we also show the ratios between the X-ray masses and 
the true masses of the clusters (asterisks) and the ratios between the 
masses determined via the hydrostatic equilibrium equation using the 
true gas density and temperature profiles and the true masses (dashed 
line). The diamonds and the triangles show the results published by 
M08 and Z09, based on the analysis of a sample of 18 and 12 galaxy 
clusters, respectively. The data-points are slightly shifted along the A- 
axis at each over-density, in order to avoid overlapping and facilitating 
the comparisons. 

the ratio M3d,xM^3d,l- Not surprisingly, the mass ratio which 
has the smaller scatter is Mid,x/M2d,l, whose r.m.s. is ~ 13%. 
However, the radial dependence of this mass ratio is stronger 
than for M3d,x/^3d,l or M2d,x/^2d,l as indicated by the vari- 
ation of the mean and of the median as a function of the over- 
density A. Obviously, the mass ratio which is affected by the 
largest scatter is M 2 d,x/M 3DiL (~ 25%). 

5. Summary and conclusions 

In this paper, we have tested several methods to derive the 
mass of galaxy clusters via lensing and X-ray observations. 
To do so, we have used two pipelines for simulating observa- 
tions of galaxy clusters obtained from N-body and hydrody- 
namical simulations. These are the code SkyLens, which pro- 
duces mock observations with a variety of telescopes in the op- 
tical bands, including lensing effects by matter along the line 
of sight, and the code XMAS, which allows to mimic X-ray ob- 
servations of galaxy clusters with the Chandra and the XMM- 
Newton telescopes. In our analysis, we have used three massive 
galaxy clusters at redshifts between z = 0.2335 and z = 0.297 '. 
Each cluster has been projected along three orthogonal axes 
and each projection has been analyzed individually. For each 
cluster projection, we have carried out both the lensing and the 
X-ray analysis using standard techniques to derive the mass 
profiles. 

Our results can be summarized as follows: 

- strong lensing parametric mass reconstructions provide 
very accurate estimates of the projected mass, but only in 



the regions probed by the strong lensing features. Our simu- 
lations show that within the Einstein rings, the input masses 
are recovered with an accuracy of < 20%. For most of the 
cases, the measured masses differ from the input ones by 
only a few percent. The precision depends on the complex- 
ity of the lens and on the number of lensing constraints 
available. For example, using central images improves the 
mass estimates in the very inner regions. On the other hand, 
extrapolating the mass models to radii larger than those 
probed by strong lensing may easily result in mass esti- 
mates wrong by up to 40 - 50%. In particular, the assump- 
tions made for modeling the star distribution in the BCG 
affect the mass profiles at large radii. Although the very 
central regions of numerically simulated clusters may be 
not well representative of those of real clusters, given the 
uncertainties on the treatment of small-scale physics of the 
intra-cluster gas, our results warn against potential biases 
which may affect the strong lensing mass estimates; 

- weak lensing methods for measuring the projected mass 
perform well with galaxy clusters characterized by regu- 
lar mass distributions, i.e. without massive substructures 
which produce a significant shear. If present, such substruc- 
tures need to be properly modeled. If not, the cluster masses 
can be over- or under-estimated by up to 50%, when assum- 
ing that the shear signal is tangential to a unique mass cen- 
ter. Non-parametric methods combining weak and strong 
lensing, which by construction do not assume any symme- 
try in the lensing signal, perform well even in the case of 
clusters disturbed by massive substructures. The mass esti- 
mates at /?2500> /?5oo, and 1?200 are typically recovered within 
a 20% accuracy, being more accurate at the smallest radii; 

- deprojection of the lensing masses introduces a signifi- 
cant scatter in the ratio between estimated and true masses. 
Cluster triaxiality causes that the accuracy of 3D-mass es- 
timates depends on the orientation of the lens with respect 
to the line of sight. For clusters whose major axis points 
towards the observer the 3D-mass is over-estimated, while 
the opposite happens for clusters oriented perpedicularly to 
the line of sight. Cluster substructures also cause errors dur- 
ing the de-projection, as the distance of the substructures 
from the cluster center along the line of sight is unknown; 

- X-ray masses are typically biased low by 5 - 20% due to 
the lack of hydrostatic equilibrium in the simulated clus- 
ters. The average bias between /?2500 and R200 is ~ 10%. In 
fact, the accretion of material from the surroundings during 
the cluster formation process causes bulk motions which 
contribute to the pressure support of the gas. As shown 
in previous works, this effect is radial dependent, i.e. it is 
larger in the outer regions. We investigated whether the ra- 
tio of X-ray-to-lensing masses a& can be used as a tracer 
of hydrostatic equilibrium. By averaging over the sample, 
we find that such ratio is indeed reproducing the lack of hy- 
drostatic equilibrium in the input simulations, although the 
number of clusters available is limited and despite the ef- 
fects of triaxiality and substructures which affect the cluster 
3D-mass estimates. 

- the gas mass is well reconstructed within the region from 
which we can extract a surface brightness profile. Almost 
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Table 4. Summary of the comparison between lensing, X-ray, and true masses. We report the mean, the r.m.s., the median, the first and the third 
quartiles, and the minimal and maximal values of the ratios between 3D and 2D masses obtained via lensing and X-ray analyses and from the 
input simulations. The lensing masses are measured with the SL+WL method. For the X-ray masses we use both the forward and the backward 
methods. 



ratio 


A 


mean 


r.m.s. 


median 


1st quart. 


3rd quart. 


min 


max 




2500 


0.9241 


0.1147 


0.8854 


0.8451 


0.9730 


0.8007 


1.1143 


M 3D™ / M 3D,t™ 


500 


0.8811 


0.0538 


0.8950 


0.8594 


0.9059 


0.7640 


0.9493 


MfgyM 3D-true 


200 


0.8947 


0.0625 


0.9162 


0.8534 


0.9293 


0.7749 


0.9759 


M 3 b g*/A*3D,lmc 


2500 


0.9205 


0.0690 


0.9018 


0.8669 


0.9669 


0.8512 


1.0407 


M 3Dx/ M 3D,t™ 


500 


0.8961 


0.0582 


0.8826 


0.8817 


0.9288 


0.7753 


0.9853 


M^/M 3D-true 


200 


0.8832 


0.0664 


0.9019 


0.8878 


0.9126 


0.7192 


0.9559 


M3D.L/M3D.lm(j 
M3D,L/*f3D,true 
^3D.L/^3D.lmc 


2500 
500 
200 


1.0233 
0.9963 
9807 


0.1686 
0.1661 
1734 


0.9683 
0.9374 
8891 


0.9093 
0.8930 
8687 


1.1033 
1.0251 
1 0783 


0.8290 
0.8362 
8234 


1.3687 
1.3438 
1 2879 


Af'°™ /M 1dl 
3D.X' JL,L 


2500 


0.9303 


0.2193 


0.8917 


0.8303 


1.0034 


0.5850 


1.3441 


3D.X' JL >> L 


500 


0.9062 


0.1521 


0.9429 


0.8837 


1.0167 


0.5686 


1.0768 


3D.X 


200 


0.9395 


0.1802 


1.0091 


0.7915 


1.0815 


0.6017 


1.1127 


M$rft/M 1DL 
3D.X' JL,L 


2500 


0.9230 


0.1776 


0.9109 


0.8416 


0.9934 


0.6219 


1.2553 


3D.X' JU . L 


500 


0.9244 


0.1731 


0.9415 


0.8926 


1.0690 


0.5770 


1.1108 


M 3Dx' M3D . L 


200 


0.9307 


0.1937 


1.0088 


0.8374 


1.0541 


0.5584 


1.1353 


M2D,L/*f3D,tme 
M2D,L/M3D.lrue 
M2D,L/^3D,tme 


2500 
500 
200 


1.4898 
1.2924 
1.2171 


0.2485 
0.2151 
0.2092 


1.4317 
1.2127 
1.1109 


1.3671 
1.1838 
1.0700 


1.5457 
1.3059 
1.3212 


1.1611 
1.0607 
1.0545 


2.0612 
1.7788 
1.6244 


MjjgJ/M 3Dmie 


2500 


1.2850 


0.0757 


1.2744 


1.2627 


1.3288 


1.1276 


1.3839 


M 2D™ /M 3D.truc 


500 


1 .0477 


0.0852 


1.0376 


0.9853 


1.0810 


0.9179 


1.1760 




200 


0.8961 


0.0835 


0.8839 


0.8446 


0.9382 


0.7449 


1.0162 




2500 


1.3211 


0.0754 


1.3317 


1.2903 


1.3520 


1.2052 


1.4415 


M 2D°X /M 3D,t™ 


500 


1.1164 


0.1121 


1.1 106 


1 .0643 


1.1270 


0.9497 


1.3130 


M^ x /iU 3Dilrlie 


200 


1.0124 


0.0990 


1.0200 


0.9919 


1.0436 


0.8108 


1.1585 


M 2D™ / M 2D,L 


2500 


0.8875 


0.1733 


0.8901 


0.8221 


0.9779 


0.5471 


1.1702 


M 2D™ / M 2D,L 


500 


0.8314 


0.1503 


0.8790 


0.7545 


0.9272 


0.5160 


0.9834 


Mf™/M 2D-L 


200 


0.7596 


0.1629 


0.7949 


0.6388 


0.8664 


0.4586 


0.9636 


M 2DX' M 2D.L 


2500 


0.9103 


0.1657 


0.9301 


0.8634 


1.0259 


0.5847 


1.1584 


M^ X /M 2D-L 


500 


0.8889 


0.1825 


0.9475 


0.8150 


1.0299 


0.5339 


1.0827 




200 


0.8609 


0.1951 


0.9259 


0.7507 


0.9638 


0.4991 


1.0827 


M 2D™ /M 3D,L 


2500 


1.2929 


0.2579 


1.3392 


1.1445 


1.5218 


0.8239 


1.6389 


M 2Dx/ M 3D.L 


500 


1.0801 


0.2064 


1.1532 


0.9612 


1.1762 


0.6830 


1.3297 


M 2D™ X / M 3D.L 


200 


0.9450 


0.2123 


0.9728 


0.7828 


1.0552 


0.5783 


1.2341 


M^*/M 3D-L 


2500 


1.3274 


0.2560 


1.3326 


1.2020 


1.5853 


0.8806 


1.6223 


M^/M 3D , L 


500 


1.1553 


0.2502 


1.2023 


1.0383 


1.3064 


0.7067 


1.4481 


M!£<*/M 3D>L 


200 


1.0711 


0.2524 


1.1663 


0.9198 


1.1748 


0.6295 


1.3760 


M 3D™' M 2D.L 


2500 


0.6407 


0.1612 


0.6212 


0.5728 


0.6556 


0.3885 


0.9597 


Mf° D ™/M 2D , L 


500 


0.6985 


0.1166 


0.7281 


0.6838 


0.7519 


0.4295 


0.8489 


M 3D r V M 2D.L 


200 


0.7557 


0.1407 


0.8147 


0.6459 


0.8509 


0.4771 


0.9013 




2500 


0.6352 


0.1302 


0.6173 


0.5895 


0.6479 


0.4129 


0.8963 


M^/M 2 D,L 


500 


0.7120 


0.1299 


0.7391 


0.7006 


0.7906 


0.4359 


0.8757 


M^ X /M 2D , L 


200 


0.7485 


0.1510 


0.8009 


0.6834 


0.8633 


0.4428 


0.8933 





indipendently of the dynamical state of the cluster the gas 
mass con be recovered, on average, with average deviation 
of 1 + 3% at /?250o and 7 ± 4% between /? 5 oo and /?2oo- 



These results indicate that the comparison between lensing 
and X-ray mass estimates in individual systems can be very 
misleading, due to the fact that these two observables probe 
different quantities. While lensing probes the matter projected 
along the line of sight to the clusters, the X-ray emission from 
the intra-cluster gas probes the three-dimensional gravitational 
potential. Thus, a direct comparison between lensing and X- 
ray masses is possible only in those systems which are char- 
acterized by large degrees of spherical symmetry. Additional 
complications are caused by substructures, which largely af- 
fect the weak lensing measurements, at least using some meth- 
ods of converting the shear into a mass estimate. Combining 
strong and weak lensing, we obtain on average a very small 
bias in the mass estimates (~ 2%) but the scatter is pretty high 
(~ 17 - 23%). On the contrary, X-ray masses are affected by 



a more substantial bias (~ 10%) but the scatter is smaller by 
almost a factor of three. 

Finally, strong lensing alone provides mass estimates which 
are reliable only in a very small region around the cluster core, 
which cannot be probed by X-ray. Extrapolating the strong 
lensing mass models to the radii probed by X-ray (and vice- 
versa) may result in significantly large mismatches between 
mass estimates. Note however that the difficulty to fit the strong 
lensing data in simulations with cooling and star formation us- 
ing the approximations which are commonly used in observa- 
tions is very likely to be due to the unrealistic core structure of 
the simulated clusters. 

Nevertheless, it is remarkable that, by averaging over the 
whole sample, we have been able to reproduce the ratio be- 
tween true and hydrostatic equilibrium masses seen in the sim- 
ulations, under the assumption that the former is given by lens- 
ing and the latter by X-ray. This suggests that, with a suffi- 
ciently large sample of clusters with wide-field imaging and 
X-ray observations, the impact of several noises affecting the 
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mass estimates can be reduced. These are positive news and 
suggest that weak lensing masses may be used as reliable cal- 
ibrators for the scaling relations involving the X-ray observ- 
ables, like the M - L x , the M - T x , and the M - Y x relations. 
Moreover, the scatter due to triaxiality in the calibration can 
potentially be reduced. In fact, although its radial dependence 
is stronger, the ratio between 3D X-ray masses and 2D lensing 
masses is affected by a smaller scatter than the ratio between 
3D or 2D X-ray and lensing masses. 
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Appendix A: Lensing definitions 

In this Sect., we summarize some lensing definitions which will be 
useful in the rest of the paper. 

We start with an isolated lens whose surface-mass density is 1.(6) 
at the angular position on the sky. Its lensing potential is 



AG A A 



) ln|0-6»'| , 



(Al) 



where A,s,is are the usual angular-diameter distances between the ob- 
server and the lens, the observer and the source, and the lens and the 



source, respectively. The reduced deflection angle experienced by a 
light ray crossing the lens plane at is the gradient of the potential, 



aid) = Vi/Kfl) 



(A.2) 



The image positions 6 for a source located at /} are given by the lens 
equation 



P = - a(9) 



(A.3) 



For sources much smaller than the typical scales on which the lens 
properties vary, the lens mapping can be linearised. The deformation 
of images with respect to the source is then given by the Jacobian 
matrix 



i= de-[ d - J -dm-j-[ -y 2 1-K + yJ- (A - 4) 



Here, k is the convergence 
!(0) 1 

Z. cr Z 

i.e. the surface-mass density scaled by its critical value 



(A.5) 



A s 



(A.6) 



4nG A A 

The distortion is described by the two components of the shear, 

Ti = ^ (^n - fe) , 72 = i/fn , (A.7) 

which are combined into the complex shear y = y x + iy 2 . We further 
use the common abbreviation 



5V(0) 



(A.8) 



The inverse of the Jacobian determinant defines the magnification 
factor 



1 



detJ?[(0) 



(A.9) 



The loci on the image plane where det^l = 0, i.e. where the magnifi- 
cation diverges, are the tangential and the radial critical lines. These 
form where the tangential and the radial eigenvalues of the Jacobian 
determinant nullifies: 



A t = l- K -y=0 
A, = 1 - K + y = . 



(A. 10) 
(A. 11) 



The inverses of the eigenvalues of the Jacobian matrix define the tan- 
gential and the radial magnifications: images forming close to the tan- 
gential or to the radial critical lines are highly elongated tangentially 
or perpendicularly to the critical lines, respectively. 

The critical lines are mapped onto the caustics on the source 
plane, via the lens equation. The caustics encompass regions of the 
source plane which are characterized by different multiplicities of the 
images. 

Outside critical curves, image ellipticities are determined by the 
complex reduced shear 



8(0) 



7(0) 
1 - k(9) 



(A. 12) 



In the weak-lensing limit, k« 1, and the reduced shear approximates 
the shear, g « 7, to first order. 

Source and image shapes are quantified by the complex ellipticity 

„ - h 

(A.13) 



a + b 
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where a and b are the semi-major and semi-minor axes of an ellipse 
fitting the object's surface-brightness distribution. The position angle 
of the ellipse's major axis is §. The expectation value of the intrinsic 
source ellipticity e s is assumed to vanish. 

A sufficiently small source with ellipticity e s is imaged to have an 
ellipticity 

- fs+g (A.14) 



1 + g'e s 



where the asterisk denotes complex conjugation. This equation illus- 
trates that the lensing distortion is determined by the reduced shear, 
which is the only lensing quantity directly accessible through mea- 
surements of galaxy ellipticities alone. 



Appendix B: Individual mass profiles 

In this Sect, we show the 3D-mass profiles (normalized to the true 
mass profiles) of each individual cluster projection in our sample. This 
allows the reader to compare the mass estimates at any radius different 
from ^2soo> and R200, which are already discussed in the paper. 
The profiles are displayed in Fig. B.l for all the lensing and the X-ray 
methods. Additionally, we also show the mass profiles derived from 
the true gas and temperature profiles under the assumption of hydro- 
static equilibrium. These last curves are quite noisy, especially in the 
case of g72, because we have not applied any smoothing to the data. 
In the case of g51 the extended cooling region in the cluster center 
causes the hydrostatic mass M H eq to be larger than the true mass. The 
X-ray derived profiles reproduce well the trends seen in the Mheq pro- 
files although some differences are present. Note that, while we mask 
many cold blobs in the simulated observations, we are not removing 
these structures from the input clusters when calculating M H eq- Thus, 
some differences between the profiles are expected. 
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Fig. B.l. Ratios between the mass profiles recovered from the lensing and the X-ray analyses and the true 3D mass profiles of each cluster. The 
meaning of each line is explained in the legend. The vertical dashed lines in each panel mark the positions of ^2500, R500, and R200, as derived 
from the lensing analysis (see text for more details). 



